I 


otf TEc $<v 


* 7 ? 

-& 5 -/ 




y o 



Rensselaer Polytechnic Institute 

Troy, New York 12181 


L 


(NASA-CB-1 36776) PRELIMINARY NUMERICAL 
ANALYSIS OF IMPROVED GAS CHROMATOGRAPH 
MODE! {Rensselaer Polytechnic Inst.) 

100 P HC $7.00 CSCL 14B 


H74-16097 


Onclas 

G3/1 4 15917 



R.P.I. Technical Report MP-3& 

PRELIMINARY NUMERICAL ANALYSIS OF 
IMPROVED GAS CHROMATOGRAPH MODEL 

Philip T. Woodrow 

National Aeronautics and Space 
Admini s tra ti on 

Grant NGL 33 -OI 8 -O 91 


Analysis and Design of a Capsule Landing System 
and Surface Vehicle Control System for 
Mars Exploration 

September, 1973 

School of Engineering 
Rensselaer Polytechnic Institute 
Troy, New York 


/ 



ABSTRACT 



In continuance of system studies on the science package for an 
unmanned Martian roving vehicle, a new mathematical model for the 
gas chromatograph has . been developed which incorporates the hereto- 
fore neglected transport mechanisms of intraparticle diffusion and 
rates of adsorption. . Moment analysis of this Inter-Ihtraparticle 
Adsorption Model has showed the model to he more capable of predict- 
ing spreading in experimental chromatograms. Because a closed- form 
analytical solution to the model does not appear realizable, techni- 
ques for the numerical solution of the model equations are being 
investigated. Criteria have been developed for using a finite term- 
inal boundary condition required in numerical solutions in place of 
an infinite boundary condition used in analytical solution techniques. 
The method of. Finite Differences appears computationally inefficient 
for application to equations of the type to be solved. The class 
of Weighted Residual methods known as Orthogonal Collocation is 
presently being investigated and appears promising. 
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I. INTRODUCTION AMD SUMMARY 

The mathematical modelling of the gas chromatorgaph is one sub- 
task of a group effort designed to define fundamental system design 
criteria necessary for an optimal design of a combination gas chroma- 
tograph - mass spectrometer which is to be part of an unmanned mission 
to Mars'. The task which must be performed by this part of a Martian 
Roving Vehicle is the analysis of samples to determine the existence of 
organic matter and living organisms on the Martian surface. The analy- 
sis will involve the subjection of gaseous, liquid, and solid samples 
to biological and chemical reactions, with subsequent product separation 
and identification using the gas chromatograph - mass spectrometer 
system. 

The chromatograph may be looked upon as a separating device where 
the phenomenon of adsorption-desorption is utilized. Owing to the 
different characteristics of various chemicals, each species will adsorb 
and desorb at different rates when exposed to a packed bed of granular 
particles with or without a liquid substrate. Because of the unique 
behavior of each chemical, a multi -component sample may be injected 
into a chromatograph and elute as separate waves of specific chemical 
species. 

The transport mechanisms which have been included in previous 
model formulations are all -interparticle mechanisms with simple adsorbed 
phase behavior assumed. These previous model formulations have proven 
incapable of adequately predicting component behavior in all cases. 
Consequently, a new model is developed which includes both interparticle 
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and intrap article transport mechanisms. This. model is analysed in the 
Laplace transform domain using the method of moments. The first three 
moments of the impulse response of the model are derived. Using actual 
input data, predictions for the first three moments of the output data 
are made and are compared with actual output data and predictions of 
a simpler interp article model. The results indicate that the new model 
is more capable of the prediction of the moments of the actual data. 

Because the mathematical complexity of the new model prohibits 
a direct, closed-form analytic expression for a response, investigation 
of numerical techniques, applicable to the equations of the old and new 
models is made. The numerical techniques require a finite terminal 
boundary condition as opposed to an infinite column boundary condition 
$ used in analytic solution (when possible) of the chromatographic model 
partial differential equations. Using a simple, transient diffusion- 
convection equation, criteria are developed wherein a finite terminal 
boundary condition can be applied to yield infinite column behavor at 
the bed outlet. • £ 

An analysis of two methods for the numerical solution of partial 
differential equations of the type encountered in the chromatograph 
modelling work is subsequently made. The technique of Finite Differ- 
ences is rejected due to excessive computer time required to produce 
model simulations. The technique of Orthogonal Collocation, while not 
established as the best method, offers promise and is the current area 
of modelling endeavor. 
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II. CHROMATOGRAPH SYSTEM MODELLING 

A. Chromatograph Modelling Background 

One area of the overall gas chromatograph systems study has 
been the mathematical modelling of the chromatograph system. Work in 
this area has been carried out by several investigators (Sliva, 1968; 
Voytus, 1969; Taylor, 1970; Keba and Woodrow, 1972). A course has been 
pursued wherein successively more complex models have been considered. 
These models have all yielded analytical expressions from which a simu- 
lated chromatogram could be computed directly. Comparison of predicted 
system behavior with actual system data has directed modelling efforts 
to consider, more adequate and hence more complicated models . 

Prior to this investigation, the most complex model pro- 
l posed for the chromatograph system was based on an interparticle phase 
mass balance and an adsorbed phase mass balance. Several transport 
mechanisms were included: axial diffusion, convection, and mass transfer 

between the inteiparticle and adsorbed phases. A linear isotherm was 
used to des< ribe the adsorption kinetics. This model has been studied 
and compared (Keba and Woodrow, 1972) for the cases of finite rates of 
mass transfer to the adsorbed phase (nonequilibrium adsorption) and 
infinitely high rates of mass transfer to the adsorbed phase (equilibrium 
adsorption). In both cases, simulations using the models failed to 
predict the degree of dispersion exhibited by many of the experimental 
data. It was concluded that additional transport mechanisms, e.g., 
intraparticle diffusion, may be contributing appreciably to the 
overall adsorption-desorption process. Hence, further model develop- 
ment and analysis was indicated. 
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B. Development of the Inter-Intraparticle Adsorption Model 

Previously, the intraparticle region of the chromatograph 

packing material has been modelled as being nonexistent or as a region 

where the transport processes occur at such a rapid rate so as not to 

significantly affect the dynamic behavior of the system. It is the 

purpose of this section to reformulate the chromatograph system model 

by including the transport processes which are presumed most likely to 

affect the dynamics of the adsorption-desorption process within the 

chromatograph packing material. 

Figure 1 presents graphically the transport processes to be 

modelled. The sample to be separated is injected into a relatively 

inert carrier gas, e.g., helium. As this slug of sample is transported 

» down the chromatograph by the carrier gas, the various species diffuse, 

adsorb, and desorb. Diffusion of the chemicals in the direction of the 

carrier gas flow in the interparticle region is represented by the 

dimensionless parameter,. Pe^, which is determined by the system fluid 

mechanics. Mass transport from the interoarticl^ region to the intra- 

particle region is represented by a dimensionless parameter, N. 

tob 

which is essentially deteimined by the system fluid mechanics. Diffusion 
in the intraparticle region is represented by dimensionless parameter, 

Pe , which is in part determined by the properties of the particle 

A 

packing. The rate of adsorption within the particle is characterized 
by the dimensionless parameter, N DTT . Adsorption- desorption within the 
particle is represented by piR^, a thermodynamic parameter peculiar to 
each species. This parameter contains an equilibrium constant, m, 



Figure 1 

chromatographic column 

INTER-INTRAPARTICLE MODEL CONCEPTS. 
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and the quantity R£. R^ is the ratio of moles of fluid within the 
particle to the moles of adsorptive sites within the particle. The 
quantity is directly related to the quantity R^ where R^ is 

the ratio of moles of fluid within the total "bed to the moles of adsorp- 
tive sites within the total "bed. The relationship between these quanti- 
ties is 



The reason for noting this relationship so that the parameter mR^ has 
been noted in previous models and the above relationship serves as a 
unifying concept for the new model formulation which follows. 

With the above concepts in mind, the following set of dim- 
i ensionless equations has been derived. based on the assumptions which 
follow: 


An inteiparticle phase mass balance: 



_ J_y 

^z 2 z 


N toG ' y “ 


= ^_y 


( 2 ) 



i 


* See Section IX, Nomenclature, 'for definition of terms. 
** See Appendix A for derivation. 
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An adsorbed phase mass balance : 

(~y = % (y i - o w 

A thermo dynamic relationship 'between the intraparticle 
and adsorbed phases : 

y i* *= mx a (5) 

The above equations are valid under the following assumptions: 

1. The column is isothermal. 

2 . The carrier gas velocity profile is flat. 

3 . • The axial diffusion coefficient is a composite factor 

which may or may not have a turbulent component. 

4. The gas composition is approximately constant in the ’ 
radial direction at a given axial position. The con- 
centration gradient occurs in a thin boundary layer at 
the interparticle-intraparticle interface. 

5. The gas composition within the particle is approximately 
constant in the angular direction at a given radial 
position; the concentration gradient occurs only in a 

■ thin boundary layer near the adsorbent surface. 

6. The adsorbent layer is so thin that there is no diffu- 
sional resistance within the layer in the direction 
normal to the surface . 

7. The diffusivity in the adsorbent layer is so small that 
there is no diffusion in the direction parallel to the 
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surface in the intraparticle radial direction. 

8. The net rate of adsorption for the carrier gas is 
negligible , 

9- Only one component is adsorbed and its gas phase com- 
position as a mole fraction is small compared to unity. 
10. The carrier gas behaves as an ideal gas. 


An applicable set of boundary and initial conditions are 


as follows: 


Initial Conditions : 
y (z, 0) .= 0 

y. (z, r, 0) = 0 

\ ( z > r ’ °) = 0 


Boundary Conditions : 
y (0, O) = 

[(f) a L (§) /PeJ |^ = 


N toG (y - y^; when r=l ... 


( 6 ) 

( 7 ) 

(8) 

( 9 ) 

(10) 


<^i 

d r 


= 0 


r = 0 


(n) 


lim y (z, 6) = finite 

Z->0D 


( 12 ) 


These conditions reflect a sample-free column at zero time, a sample 
injected as an impulse, mass transfer between the interparticle and 
intraparticle regions, no concentration gradient at the center of the 
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column packing, and no end effects at the column exit. 

For the systems under consideration it has been shown by 
Keba and Woodrow (1972) that inclusion of the parameter N^ oG is of 
minor importance. If one were to consider the case of infinite rates 
of mass transfer, i.e., N^ — ►qd , the coupling condition given by 
equation (10) would be replaced by 


y i (z, 1, O) = y (z, 6) 


(13) 


Ulus, a. model in the form of a set of coupled, partial 
differential equations is proposed. Prior to consideration of the time 
domain solution of the equations, a moment analysis can be made to 
ascertain the predictive capabilities of the proposed model. This 
analysis is the subject of the next part of this report. 


III. MOMENT ANALYSIS OF THE INTER- INTRAPARTICLE ADSORPTION MODEL 


A. Theory and Background 

An analysis of a proposed model can be made prior to deter- 
mination of the model's time-domain solution to yield the gross charac- 
teristics of the impulse response of the model. In addition, because 
of the poor predictions of previous models with respect to chromatogram 
spreading, it is desirable to know the nature of the response of the 
proposed model for the pulse-type forcing functions used in experi- 
mental work. The nature of the response can be characterized by statis- 
tical quantities known as moments which may be obtained without knowledge 
of the time- domain model solution. The moments may be derived directly 
from the Laplace domain solution of the model. The following develop- 
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ment will indicate how the moments of a model are obtained and how 
the analysis can be extended to give the moments of systems forced by 
general pulse -type inputs. 

The impulse response of the chromatograph may be viewed 
as the residence time frequency distribution (Douglas, 1972). This 
quantity resembles the probability distribution function which appears 
in statistical analysis. The moments of the distribution function 
about the time origin are defined by the following: 


'00 


-00 


= 


O n f (O) d -9 / 1 f (O) d 9 


(14) 


'o 


'0 


where 

f (O) = the distribution function being analysed. 

The denominator of equation (l4) is the area under the function. The 
relationship of the moments about the origin to the Laplace transform 
is developed in Appendix B. The result is: 

n 


d. 


= (-l) n limi ^ — [f (s)l ( / lim f (s) 

n ~ L J ‘ s-j>0 


(15) 


s->0J ^ n 
O s 


where 


7 (b) = l [f (e)] 


e' s9 f(0) d 9 


(16) 


Interest also centers on the moments about the first abso- 
lute moment or mean, Mathematically these moments are defined by: 

r „ r 00 

(17) 


= 


o-' 


n r s 

(e - P x ) f («) d e/ f ($) d O ; n - 


O 



11 . 

These moments about the mean are directly related to the moments 

about the origin. The relationships are obtained by formal expansion 
of equation (17). Appendix B gives the relationships for n=2 and 
n=3. For n=2, the moment about the mean is exactly the variance of 
the response. For n=3> the moment about the mean is related to the 
skew of the response. 

One can use the preceding to develop equations relating 
the moments of system responses for arbitrary pulse-type forcing 
functions (see Appendix B for details). That is, given the system 
input data (the moments of which can be computed from equation (lU)) 
and the system transfer function (the Laplace transform of the impulse 
response), the moments of the system response may be determined and 
V compared with the moments of the actual output data. Referring to 
the block diagram in Figure 2, the results are: 


= 

V A c 


( 18 ) 

^1Y = 

"is + 

^1G 

(19) 

"2Y = 

"2X + 

^2G 

(20) 

"3Y 

^3X + 

m 3G 

(21) 

Equation (18) states 

that the 

area under the output curve 

is the product 


of the area under the input curve and the impulse response curve. 
Equation (19 ) states that the mean of the output occurs at the sum of 
the mean of the input function and impulse response. Equation (20) 
states that the variance of the output is the sum of the variance of 
the input function and the variance of the impulse response. Equation 



12 . 



Y(s) “ L £y (e)^| ; y (0) is the system response 
G(s) = system transfer function 


e 


Figure 2 . Typical System Block Diagram 
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(2l) states that the third moment about the mean of the output is the 
sum of the third moments about the means of the input function and 
impulse response, respectively. 

This technique can also be used for estimating ^rtem 
parameters. Douglas (1972) uses an equation similar to equation (20) 
to estimate an axial Peclet number for a packed bed. Schneider and 
Smith (1968) apply moment analysis to estimate adsorption equilibrium 
constants, rate constants, and intraparticle diffusivities for a 
chromatographic system modelled similarly to that. of Part II. . 
However, accurate parameter estimation using this method is limited 
by the accuracy of the data used for analysis. 

B. Application of Moment Analysis to the Inter- Intraparticle 

Adsorption Model 

The previous section outlined a method which can be used 
to analyse pulsed systems to determine the gross characteristics of 
the system response. This section will document an application of the 
concepts of moment analysis to the proposed model of Part II . ’ 

Consider the set of partial differential equations, boundary 
conditions, and initial conditions, equations (2) through (12). A 
Laplace transform domain solution for the impulse response or transfer 
function was derived and appears in Figure 3j details appear in 
Appendix C. 

Applying the definition given by equation (15) and using 
equation (17), the moments and are derived for the 

impulse response of the Inter- Intraparticle Adsorption Model. The 





Figure 3 


where : 


Y(1>b) 

7(») 

V(s) 


a l( s ) 


exp 




Pe 


E 


7(b) pe 


E 


N t0 , (1 - X (s) ) + s 


*b sinh 


inh (^) 


Qh-I )sinh (-^a”) +\fs^ cosh (^a^)l 


/•*- 

* J. ' J. 



-hL mR_ 

2 

- 

RU I 

<T) Pe A 

s.' 

+ W RTT + S 

( st Vi) . _ 


N 


toG 


3 P (1-e) /L% _1_ 

e V R ; Pe. 

A 


a 


Particle porosity 
Bed void fraction 


Transfer Function for the Inter- Intraparticle Adsoiption Model. 
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results are presented in Figure 4; details of the manipulations appear 
in Appendix D. 

The parameters Pe , N, _, and Pe. can he estimated a 

■ill UOG- A — 

priori ♦ The parameters mR and N are not predictable a priori . 


Previous modeling analysis has estimated mR^ by a curve fitting 

process (Benoit, 1971 ). The estimation of IL will most likely 

nU 

involve curve fitting also . 

An analysis was made using existing single component data. 

The parameters Pe_, Pe., and N, were estimated using existing 

Ji A 0O(j 

correlations. The values of mR Q which were estimated by Keba and 

Woodrow (1972) using simpler models were used and the parameter N 

was varied. Tables 1 and 2 give results of this analysis for acetone 

at. 100° C and ethylene at 50° C. Both experiments used Chromasorb 102 

column packing a porous material. In each case, the moments for the 

impulse response of the model were computed using the equations given 

in Figure U. Use of system input data and equations ( 19 ) through (21) 

give predictions as a function of N DTT f^r the output moments. These 

i\U v 

predicted values are compared with actual moments of the output data 

and with the predictions of the simpler, interparticle equilibrium 

adsorption model. Expressions for the moments of the simpler model 

were initially developed by Voytus ( 1969 ). 

The results indicate that the proposed model can more 

closely predict the characteristics of the output data than the simpler, 

interparticle model. The results indicate that a value of K on 

RU 

the order of several hundred will give a predicted second moment very 





Hi = 1 + + ^ 1 " € ^ ^/ € 


= 2 


(^) 2 /p< 


'E 


+ 2 


Ql-e) p/e"] |(l + 1/mRj) 2 [(b/l) 2 Pe A /l5 + (1-e) p/e tr toG ]+ l/W^ 


(mRj)' 


^3 = 6 H 2 /Pe E + 6 ["(1-e) p/e[] jj 5 Q 1 + (mB r ) 2 3 [(b/L) 2 Pe A /l5 

•■■■.' + (i-e) P/e N toG ] + (1 + 1/mR!) 3 [((1-e) f/e sj 2 

+ 2 (1-e) p (r/l) 2 Pe A /l5 e - 23 (R/t) U Pe A 2 / 3X5~J' 

+ 1/H^ 2 (mR-j.) 3 } 


Figure 4. Moments of the Impulse Response of the Inter- Intraparticle Adsorption Model. 
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_^0_ 

O.029 


TABLE 1 

MOMENT ANALYSIS AND PARAMETRIC 
STUDY - ACETONE 100 °C . 




observed 


„ M- 

predicted 


^L, 


( 2 ) 

predicted 


173.29 


158.69 


156.49 




observed 


_ (D 

^2, predicted 


( 2 ) 

^ 2 , predicted 


100 

200 

400 

800 

1600 

3200 

64oo 

12800 

25600 


815.67 


V 



4 


3 , observed 


100. 254 o 4.0 

200 
4oo 
800 
1600 
3200 
64oo 

12800 w. 

25600 -JL 


977.55 

723.41 

686.34 

517.80 

W3 . 53 

466 . 1*0 

457 . 83 . 

453-55 

451.41 


437.28 


V 


4 


(1) 

3, predicted 

23192.2 

20454.8 

19745.3 
19555 . >■ 
19501.7 

19485.1 

19480.4 

19477.2 

19476.2 


. ( 2 ) 

l 3, predicted 

19499.2 


V 


Pe. • = 8689 . 

■E 

N too = 88960 

(L/R) 2 /Pe A = 328.2 

(1) Inter- Intraparticle Adsorption Model 

( 2 ) Interparticle Equilibrium Adsorption Model 
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TABLE 2 


MOMENT ANALYSIS AND PARAMETRIC 
' STUDY - ETHYLENE 50°C . 





'1, observed 



(i) 

predicted 



( 2 ) 

predicted 


0.194 


26 .475 


25.986 


23.719 



^2 

observed 


-d) 


“predicted 


-( 2 ) 

predicted 


;ioo 

.:200 

4oo 

800 

1600 

3200 

64oo 

12800 

35600 


7.024 


v 


13.283 

6.973 

3.817 

2.2h0 

1.451 

1.056 

0.859 

O.76O 

0.711 


0.388 


V 



^3 

^observed 


-d) 


^predicted 


"( 2 ) 


^predicted 


100 

200 

400 

800 

1600 

3200 

6400 

12800 

25600 


19.623 


V 


13.049 

3.519 

1.058 

0.403 

O.219 

0.163 

0.144 

0.137 

0.134 


0,191 




Pe E = 97^4 

"tOG = 

(1/R) (l) 2 /Pe. = 436.2 

A 

(l) . Inter- Intraparticle Adsorption Model 

(2-) Interparticle Equilibrium. Adsorption Model 
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close to the second moment of the output data. This magnitude of 

N . is consistent -with the values of N which were obtained in 
KU RU 

independent research by Schneider and Smith (1968). Tables 1 and 2 
further indicate that matching of the third moments would give different 
values , of fiL TT . However, the use of third moments is not as reliable 
because data inaccuracies are further magnified in the analysis. 

It should be noted that if one accepts the value of N_ TT 
as being on the order of several hundred for each case, all other 
parameters, excluding mR^, are of the same magnitude. The key to the 
difference in the two component behaviors is the parameter mR^. 

IV, TERMINAL BOUNDARY CONDITION ANALYSIS 

Mathematical modelling of chromatographic systems commonly 
l require solutions to equations of the form: 

(l/Pe) fe 2 y/c> z 2 ) -dy/dz - R A =^y /^0 ( 22 ) 

Application of analytical, techniques to the above equation, when possible, 
commonly utilize the terminal boundary co ditioni. 

13 m y (z, ■©) = finite ; •© > 0 (23) 

Z“*-CD 

Use of the above boundary condition in analytical work yields a great 
deal of mathematical simplification. In addition, the use of this 
boundary condition is consistent with the theory which has been devel- 
oped for prediction of the dispersion in packed beds; see, for example, 
Gunn ( 1969 ) . 
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However, when numerical techniques must "be applied to solve 
equation ( 22 ), the terminal “boundary condition given “by equation ( 23 ) 
must be replaced by a terminal boundary condition which is both computa- 
tionally expedient .and physically meaningful. A finite terminal boun- 
dary condition which has found general usage in chemical reaction 
engineering problems (Wehner and Wilhelm, 1956 ) is: 

d y ( 1 , «) /<) Z = 0 ; « > 0 (24) 

Bastian and Lapidus (195&) considered the case where R in equation 
(22) was an adsorption term. A linear relationship was assumed to 
describe the adsorption kinetics. For a step-input and the conditions 
chosen, Bastian and Lapidus showed that finite column calculations, using 
l equation (24) as a terminal boundary condition, closely approximated 
infinite column calculations, using equation ( 23 ) as a terminal boundary 
condition. , 

The analysis of chromatograph systems for pulse-type forcing 
functions has prompted consideration of t‘ e two ierminal boundary con- 
ditions. The question arises as to how the use of a finite terminal 
boundary condition -affects output prediction as compared to the infinite 
column case when the system is forced by pulse-type functions. It is 
desirable for the two predictions of column outlet behavior (z=l) to be 
similar so that the use of a priori estimates of Pe are valid in 
complicated models having the form of equation ( 22 ). 

In order to answer the above question and to establish the 
conditions under which a finite terminal boundary condition can be used 
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to yield infinite .column behavior at the column outlet (z=l), two 


relatively simple problems can be considered: 


Case I: 


(l/Pe) - R a = ° 

(25) 

y (z, 0) = 0 ; z > 0 

( 26 ) 

y (o, «). = & (©) ; © ? o 

(27) 

lim y (z, 6) ■ « finite 
z-4oo 

( 28 ) 

ii 

o 

(29) 

and 


Case II: 


(l/Pe) (a 2 y/d z) - (dy/dz) - E a = 0 

(30) 

y (z, 0) = 0 ; z > 0 

(31) 

' y (o, ©) = cf (©-) ; © ^ 6 
Oy (z Q > «) /dz = 0 ; © > 0 , z 0 -1 

(32) 

and arbitrary 

(33) 

ii 

o 

(34) 


Case I considers the unit impulse response of the simple, one- 
dimensional, axial dispersion- convection model in an infinite column. 

Case II considers the unit impulse response of the simple, one -dimensional, 
axial dispersion-convection model with the finite column boundary condi- 
tion. It is desirable to determine the conditions under which the two 
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responses are equivalent. These conditions can be determined without 

resorting to the comparisons of the analytical solutions for each case, 

through use of the method of moments. 

At a dimensionless length of unity, the column outlet, the 

* 

Laplace transforms of the two solutions are : 

• Case I t ' 

y (1, s) = exp [(Pe/2 ) - (arg)J (35) 


Case II: 


y(l, s) - exp 


(if) j(lf + arg ) ex P • arg ] - 

exp jjl-z 0 ) argjj> |(^ + argj exp (arg)J 

jj z 0 (arg)|J ... 


~ (~2 " ar S| exp 


•where 


arg = 


p ** * 

\JPe /k + Pe s 



(36) 

(37) 


Each, respective output curve can be characterized by its moments . 
Two moments are considered here - the first moment about the origin and 
the second moment about the mean. The first moment about the origin 
gives the time of appearance of the mean of the output curve. The 
second moment about the mean gives the variance of the output curve. 

These moments, as has been previously noted in Part III, are directly 
obtainable from the Laplace transform domain solution. The general 
relationships were given in equations (l4) through (17 )• Using these 


* see Appendix E for details. 
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relationships, the Case I and Case II transfer functions were analysed 
to yield: 



^2 


II 

II 


= 1 

a 2 / Pe 

= 1 + j~exp (-Pe Zq) - exp (Pe - Pe » Zq)J /Pe 

= 2/pe + exp (Pe - Pe • Z Q ) f~4/pe - 4 Z^Pe 

- 4/Pe 2 J + exp (-2 Z Q Pe)/Pe 2 

- exp (2 Pe - 2 Z Q Pe)/Pe 2 ... 


(38) 

( 39 ) 

m 


( 41 ) 


♦ If one considers the limit of the Case II moments as Z A becomes very 

i -* ' . 

large, the two results are equivalent, or: 

lim pu = l-Li = 1 

Z^od X II .1 • 

, • C 

and 

lim ]I 2 . = ~ 2 = 2/pe 

Z 0 -^od II I 

Table 3 summarizes the results of parameteric studies of the 
two moments considered for each case. The errors in Case II versus 

• Case I moments for Z Q = 1 are significant for low Peclet number. 

The error diminishes with increasing Peclet number. This confirms the 
qualitative conclusions of Friedly (1972) for high values of Pe. 



Table 3. Case I and Case* 


Errors at Zq « 1.0 


Pe 

Absolute 

Error 

W i^i n } 

Relative Error, $ 

/ ^1 I~^l II 

x 100 

\ ^ I / 

Absolute 

Error 

(m>2 II ^ 

2 

0.4323 

1+3. 23 

1.245 

4 

0.2454 

24.54 

0.3125 

8 

0.1250 

12.50 

0.07813 

16 

0.06250 

6.250 

0.01953 

; 32 

0.03125 

3.125 

0.004883 

64 

O.OI563 

1.563 

0.001221 

128 

O.OO78I 

O.78I 

0.0003052 

256 

0.00391 

0.391 

O.OOOC763 

512 

0.00195 

0.195 

0.0000191 

1024 

O.OOO98 

0.098 

0.0000047 

2048 

0.00049 

0.049 

0.0000012 

4096 

0.00024 

0.024 

0.0000003 

8192 

0.00012 

0.012 

0.0000000 


* 

^1 1 

- *1 

11 

< 

CO 

1 

0 

H 

-** 

U 2 I 

H 2 

11 

< 

10- 8 


Comparison Results 


Relative Error, $ 


l"* A 2 III 

11 x 100 

Mg 1 / 

Safe Z Q 

(Ul !!•)* 

Safe Z Q 
("^2 1^2 II 

124.5 

9.791 

II.768 

62.9 

5.254 

6.021 

31-2 

3.043 

3.328 

15-6 

1.978 

2.073 

7.91 

1.467 

1.490 

3.91 

1.223 

1.222 

1.99 

1.106 

1.044 

0.976 

■ 1.050 

1.002 

0.489 

1.024 

1.0005 

0.241 

1.011 

1.0001 

0.123 

1.005 

1.00002 

0.0615 

1.002 

1.000005 

0.0 

1.001 

1.000001 
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Table 3 also gives the value of which, when used in Case II, will 

yield output characteristics the same as Case I output characteristics. 
This means that for a given Peclet number, application of equation (33) 
at the noted Zq, will yield output characteristics at Z=1 that are, 
for all intents and purposes, the same as those predicted by Case I. 

Table 4 presents some typical values of the Peclet parameter 
for several systems. For chromatographic systems, the range of the 
Peclet number is on the order of 5,000 to 10,000. Thus in this research, 
it appears that use of the zero -derivative condition (equation 33) at 
the column exit will not cause serious problems. 

In conclusion, the comparison of the mean and variance for 
impulse responses at Z=1 for the two different boundary conditions 
has yielded guidelines which are useful when approximating infinite 
column behavior using a finite terminal boundary condition. The use of 
the criteria for- general pulse-type forcing functions would yield 
results wherein the absolute errors between the two cases would be the 
same but the relative errors between cases would decrease. The appli- 
cation of the results for models including other transport mechanisms 
(R ^ 0) may be somewhat conservative. $hen applicable to more compli- 

-A. 

cated models, the method of analysis used here will give more definite 
guidelines for each specific situation. 

V. EVALUATION OF NUMERICAL TECHNIQUES - FINITE DIFFERENCE 

This section and the one that follows present evaluations of two 
techniques which are available for the numerical solutions of the type 
of partial differential equation models that have been used and that 



Table 4. Peclet Numbers for Four Typical systems 

System Pe 

Micro Gas Chromatograph Column 233 

(Water in Helium) 

Typical Gas Chromatograph Column 5^22 

(Water in Helium) 

Typical Gas Dehydrator 1777 

■ (Water in Helium) 

Small Experimental Reactor 
(S0 2 in Air) 


155 
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are /being postulated for the chromatographic system. The two methods 
considered are the Finite Difference method and the subclass of 
Weighted Residual methods known as Orthogonal Collocation. 

Finite Difference Method 

Finite difference approximations have predominantly been used in 
the analysis of partial differential equations. To obtain numerical 
solutions to partial differential equations, one replaces the contin- 
uous variables with discrete variables. The relations between these 
discrete variables in the method of. finite differences are called finite 
difference equations. The relationships are based on Taylor series 
representations of the dependent variable. The domains of the independ- 
ent variables that are discretized form a system of grid points. Figure 
5 shows a grid representation for the transient analysis of a system 
with one spatial independent variable. The spatial dimension, Z, is 
shown as being bounded and the time variable, •©, is shown with no 
particular bound. The grid is a fixed gr*d; i.e^, spatial discretiza*^ 
tions and time discretizations are uniform for each domain. Note that 
the value of Z, the continuous space dimension is given by: 

Z = i . ( A Z) 

where i refers to a particular spatial grid point and AZ is the 
spacing between spatial grid points. Similarly, the value of •©, the 
continuous time variable is given by 

$ = 


3 • ( A «) 




29 . 


where j refers to a particular time grid point and A G is the 
interval between time grid points. 

For parabolic problems (as is the case for the second-order 
chromatograph system models), the two-level implicit method known as 
the Crank -Ni cols on method is probably most popular and is well docu- 
mented (Lapidus, 1962). In this method, the following approximations 
are made for the first and second spatial derivatives and the first 
time derivative: 


Qy/az).. ^ 1/2 


y i+i,j’ y i=i,d 
2 (A Z) 


y i+l ” y i-l, j+1 
2 (AZ) 


(dy/iZ 2 ).. » 1/2 


y i+l , i " 2sr i,j' + . y l-l,3 y i+l Jl 1+l“ 2y l,o+l +y l-l,i+l 


(A Z )' 


(A Z )' 


(a y /ae) 1;J (y ii3+ i - y ld ) / A 6 

where the i subscript denotes a coordinate in the spatial domain and 
the j subscript denotes a coordinate in the time domain. 

Preliminary studies have been made applying the Crank-Nicolson 
method to the problem: 

(l/pe ) (<J 2 y/<)z 2 ) - £y/dz = 3y/Js 

y (z,o) = 0 ; z - 0 

y ( 0 , ©) = 0 (©); © > 0 

d y (z 0 , ©)/d © = 0 ; © > 


0 
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Simulations were made with the following conditions . 

1. 0 (#) was a triangular- type pulse of duration 0.01 and 
with unit area. This is quite a sharp pulse as far as 
typical chromatograph input pulses are concerned, hut 
it was used mainly in the interest of saving computer 
time . 

2. The Peclet number was fixed at 8,000. 

3- The time increment, AO, was held at 0.0004. 

4. The response was studied at Z = 0.05. This is a drastic 
reduction in the normal spatial coordinate studied, hut, 
again, this was done in the interest of conserving 
computer time. 

5. The terminal boundary condition was applied at Z Q - 0.20. 

6. The spatial increment, AZ, was varied in the following 
sequence: 

0.0002, 0.0004, 0.0010, 0.0025 

For spatial increment values of 0.0010 and less, the simulations were 
stable. However, when A Z was increased to 0.0025, instability in 
the form of oscillation in the response was exhibited. The very small 
A Z required is directly attributed to the Pe value used. This 
instable A z value is not quite as small as the value that is pre- 
dicted by the stability criteria of Price, et. al. (1966). 

The simulation for. spatial increments of 0.0002, 0.0004, and 
0,0010 gave reasonable results when compared to results convolving 
0 (O) with the analytical impulse response. The discrepancy between 
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"between the analytic and numerical computations appeared in the magni- 
tudes of each response point - the numerical results were on the order 
of 20$ too low. This in turn affected the area beneath the response 
curve for the numerical results - all areas were on the order of 0.80 
as compared with the correct area of 1.0. The area under the analytical 
response curve was O .96 which is tolerable considering the sharp input. 
This discrepancy in response area can be resolved by adding additional 
parameters to the difference equations to yield an exact conservative 
relationship (Rogers, 1973): 

~ i ■ “■ 

System Input - System Output over 
' the interval 3 to j + 1 

N N 



i=l i=l 


where N is the total number of spatial points. This analysis was not 
performed because it was felt that the method already suffered from a 
more alarming feature - the high degree of spatial discretization which 
was necessary for the large Pe values encountered in chromatographic 
systems analysis. Extrapolation of the computing time required for the 
simulations performed yields an estimate of one to two hours of com- 
puter time required for complete simulations over the space interval 
(0, 1,0 + ). The time would naturally increase when broader input pulses 
are used. Similar conclusions on the use of finite difference schemes 
were . reported earlier (Pfeiffer, 1972). 
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Because of the high degree of spatial discretization required 
by the finite difference method and the subsequent high cost of 
computer simulations, it is felt that further pursuit of finite 
difference formulations for problems similar to the above is not 
warranted at this time and that efforts must be directed to other 
methods . 

VI. EVALUATION OF NUMERICAL TECHNIQUES ORTHOGONAL COLLOCATION 

A recent text (Finlays on, 1972) has dealt with several approxi- 
mation techniques for the solution of the differential equations which 
arise in the analysis of transport phenomena. A group of approximation 
techniques has been designated the Method of Weighted Residuals (MWR). 

A subclass of MWR is the Method of Orthogonal Collocation. This method 
« has been successfully applied to several problems in the realm of 
chemical reaction engineering. Investigators in this area include 
Ferguson and Finlayson (1970), Finlayson (1971); Villadsen and Stewart 
(1967), Villadsen and Sorensen ( 1969 ), and Villadsen (1970). The 
purpose of -"iis section is to present a summary of the theory behind 
the method and to discuss investigations that have been made relative 
to its applicability to the types of problems that must be solved in 
conjunction with the modelling of the chromatographic system. 

A. Theory and Background 

The Method of Weighted Residuals approach to the solution 
of partial differential equations starts with a representation of the 
dependent variable, y, by a finite sum of trial functions An 

example might be: 
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N 

y (z,e) = 0 O (z,«). + ^ ^ («) Pi.x (z) ( 42 ) 

1=1 

where 0 Q (Z,$) is a function which may he chosen to satisfy one or more 
boundary conditions. The functions P^(z) are normally specified and 
the time-varying coefficients, a^(o), are determined in a manner to 
give the "best" solution of the differential equation. 

The next step in the MWR is to manipulate the differential 
equation such that one side, say the right hand side, of the equation 
is zero. Then, the trial function expansion is substituted into the 
left hand side. This substitution of the trial function expansion into 
the manipulated differential equation forms what is termed the residual, 
Res, If the trial function were exact, the residual would be zero. 

In MWR, the coefficients, a^-0) are determined by specifying weighted 
integrals of the residual to be zero; 

J W.. (Res) dV = 0 ; j = 1, 2, ... I (43) 

V 

The choice of weighting functions, W. > determines what 

J 

class of MWR is to be applied. For the general collocation method, the 
weighting functions are chosen as displaced Dirac delta functions: 

w. = £ (z - z ) ; j = 1| 2, ... N W 

tJ J 

Substitution of equation (44) into equation (43) gives the result of 
forcing the residual to be zero at N specified collocation points. 

As the degree of approximation is increased, the residual will be forced 
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to be zero at an increasing number of points .in the spatial domain and 
the trial function should converge to the true solution within a given 
accuracy. 

Within the class of collocation methods is the subclass of 
orthogonal collocation. The distinguishing feature of this method is 
that the trial functions, P^(z), are chosen as orthogonal polynomials 
defined by the following relationship: 

W(Z) P. (Z) P. (z) d z * C, S a . (45) 

l 3 i iO 

a 

where [a, bj is the interval of orthogonality, w(z) is a positive 
weighting function on [a, b] , Ch is a scale factor, and <$\.. is 
the K-ronecker delta. The group of polynomials defined by equation ( 44 ) 
<±s said to be orthogonal on the interval £a, b] with respect to the 
weighting function W(Z). 

The N collocation points are chosen as roots to F^(Z), 

which is the polynomial of the next highest order in the trial function 

C ' 

expansion, the highest being P^. ^ in equation ( 42 ), The basis for 
choosing the roots of the polynomial as the collocation points instead 
of equidistant points in the interval of interest can be found in the 
theory of polynomial interpolation. Several results, as documented 
by Lanczos (1956) are summarized here: 

1 . Polynomial expansions are justified due to the funa- 
mental theorem proved by Weierstrass in 1885 which 
establishes that any continuous function in a finite 
interval can always be approximated to any degree of 
accuracy by finite power series. 
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2. The Weierstrass theorem does not imply that an 
approximating polynomial can he obtained by using 
equidistant points. This behavior was studied by 
Runge in 1901 who showed that equidistant inter- 
polation of some very simple analytical functions 
could in certain regions yield very erroneous 
results which did not disappear with increased points . 
This behavior is termed the ,T Runge phenomenon." 

3 , The difficulties which occur with equidistant inter- 
polation disappear when the zeros of the first 
neglected polynomial in the polynomial approxima- 
tion are used as interpolation points. However, 
this introduces the need to know the roots of the 
particular polynomial. 

B. Problem Formulation Using Orthogonal Collocation 

The solution of parabolic partial differential equations 
using orthogonal collocation requires several steps which are independ- 
ent of the particular equation under consideration. This section 
presents two formulations which are theoretically equivalent but which 
differ in computational and coding advantages. The first formulation, 
although somewhat more complex from a coding point of view, will be 
shown to be superior for computations. 

A trial function has been proposed, Finlayson (1972, p. 105), 
for second order systems on the spatial interval [ 0 , 1 
analysis, the trial function is of the form: 


For transient 
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N 

y (z,o) = f (e) + g (e) z + z (i-z) V a. ($) p. (z) (U 6 ) 

\ C ^ 

i=l 

The above equation has N+2 unknowns: the functions f(£), g(6) and 
(-0) , i=l, N. These are determined by the boundary conditions at 

Z=0 and Z=1 and by performing collocation at the N roots of P^(z). 
Thus, one has a set of N+2 points: 


z i - 0 

Z n+2 = 1 


and 



0=2, N+l; the roots of P^(z). 


Now, if one were to construct the approximate solution at 


these N+2 points, a matrix equation would result: 


y(z x ^) 


1 Z 1 Z i( 1 “ Z i) p 0 ^ Z i) - * * Z i^ 1 “ Z i) p W -i( Z i) 

f (e) 

y(z 2 ,e) 


1 Zg ^(l-Zg )Pq(2v-,) ’ * * Z 2^~ Z 2^ P N-1^ Z 2^ 

g (©) 

• 


• • • . 

a x («) 



1 ViW^^iVVi^'^+iH+iV^w^ 

• 

y(Z H+2' e; 


1 Z N+2 Z N+2^ 1 " Z N+2^ P 0^ Z N+2^ * * * Z N+2^ 1 “ Z N+2' )P N-l^ Z N+2 ^ 



Now define the following quantities: 


A 

y ~ 


yCz^e) 

y(z 2 ,©) 

« 

y ^ Z H+l' 6 ^ 


( 48 ) 
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A 


R 


1 2 i z 1 (i-z 1 )p 0 (z 1 ) 

1 Z 2 z 2 (i-z 2 )P 0 (.z 2 ) 


^+1 

^N+2 Z N+2 < ' 1 " Z N+2^ :P 0 ( ' Z N+2— 


VHVA 1 

z 2 (i-z 2 )p n _i ( z 2 ) 


^+1 ^ 1_ ^+1 ^ P N-1 ^ ^+1 ^ 

Z N+2 ( i-Z^g ) p H _i ( Z N+ 2 ) 


(h9) 


f (e) 

g («) 

a^e) 


^(0) 


(50) 


Use of equations (48), (49) and (50) reduces equation (47) to the more 
compact form: 

y = R f (51) 


The spatial derivatives may he expressed in a similar form: 

U -- M 


(52) 


. 2 

| = R.2 ‘ f 

(3Z — 


( 53 ) 
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where 



i = 3> N+2 



and 


R2J = [R2J o,2 = 0 ; o - 1, N+2 

flffi") d,i=z ( 1 -Z )^ 2 p i- 3 (Z j ) + 2 ( 1 - 2 Z P i- 3 (Z d 

W -J J J Q 0 


. ^ Z 

- 2 Pj.3 (Z^.); a = 1, K+2 
i = 3, If+2 


d z 


(57) 


The time-varying vector f may "be eliminated from equations (52) and 
(53) by premultiplying equation (51) by the inverse of R, if 1 , or: 

f = g" 1 y 

and 


(58) 

(59) 

c 

Equations (58) and (59) thus- yield expressions for the first and second 
spatial derivatives at the N+2 points in terms of the solution at ' 
the N+2 points . 

Alternative to the formulation of above is a formulation 

which is presented hy Finlayson. (1972, pp. 105-106). Expansion 

of equation (K6) yields an (N+l) order polynomial: 

N+l 

y (Z,0) = f (S) + £ d i (3) Z 1 

i=l 


£ 2 

= R1 R y 

d 2 ~ “ 


2 

<5 I - 

P = R2 R ' y 

d * = = “ 


(60) 
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Writing the approximate solution at the N+2 points yields a matrix 
equation similar to equation ( 51 ): 


y = § d 

where; 


1 

Z 1 

h 2 

N+l 

• • i ^1 

1 

• 

Z 2 

• 

£ 

V 

• 

N+l 

■ • • ^2 
• 

• 

1 

• 

^+1 

4i 

• 

N+l 

* * * n+l 

1 

^N+2 

Z N+2 

7 n+i 

* ’ - N+2 


£ 


(61) 


( 62 ) 



(63) 


The first and second spatial derivative vectors can he written as: 


yf = £ * (64) 

and 

^ ( 6 5 ) 
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where 



(i-1) z/- 2 
(i-D (i-2) 


' j=l, N+2 
5 i-1, N+2 


. i-3 . 3=1, 

3 . ’ 1 = 1 , 


N+2 

N+2 


( 66 ) 


(67) 


As In the first formulation, the time- varying vector, d, may he elim- 
inated from equations (6U) and (6 5) hy p re -multiplying equation (6l) 
hy the inverse of ‘ Q~‘ L , or: 

d = 3 1 y 

and 

c) y -1 

■ p = ^ £ 2 (68) 



S| S ’ 1 y 


( 69 ) 


Thus, equati ns (68) and ( 69 ) give expressions viiich are identical to 
equations ( 58 ) and (59)» Th e matrix product El R ^ is equivalent 
to and Rfr B"* 1 " I s equivalent to CJ2 Since the 

computations of Q, , ££L, and _§,2 only require knowledge of the colloca^ 
tion points and not the knowledge of the particular polynomial coef- 
ficients "being considered, one might conceivably prefer the second 
fomulation. Both formulations require the computation of the inverse 
of an (N+2) square matrix. 
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Computationally, it is desirable . for the matrix "being 
inverted to be well- conditioned with respect to inversion. An analysis 
has been made comparing the inversion qualities of the matrices R and 
£* The ease of inversion is measured by the condition number of R and 
£ respectively and with the number of decimal digits which are left 
unchanged following interative improvement of the initial Gauss -Jordan 
reduction of each matrix. Stewart (1973) discusses the problem of 
ill - conditioning and the use of iterative improvement in matrix inversion. 
Table 5 compares the inversion characteristics of R and £ for 
increasing N. The condition numbers cited are lower bounds on the true 
condition numbers relative to the noun. Appendix F shows how the 
lower bound and the upper bound on the condition number is computed. 

Except in the analysis of £ for (N+2) ^ 22, there were no practical 
differences in the lower and upper bounds. 

Table 5 indicates that the matrix R is well- conditioned 
with respect to inversion using the double-precision word length avail- 
able on the IBM 360/50 computer. In all eases, £he computation of ’ 
the product R R 1 yielded a matrix whose off-diagonal elements were 
less than or equal ,to 10 The table also shows the progressively 

poorer conditioning of £ with respect to inversion. The (26 x 26) 
case is so ill-conditioned that inversion using the available computer • 

* The noim of an (n x n) matrix A is defined as; 
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Table 5» Comparison of Conditioning of R and ^ Matrices with 
Respect' to Inversion 


Lower Bound # Lower Bound # 

. Matrix Size- on Condition of R IDGT on Condition of Q IDGT 

- K . - <4 


( 3 

X 

3 ) 

0.120 

X 

10 2 

15 

0.240 

X 

10 2 

15 

( 4 

X 

4 ) 

0.328 

X 

io 2 

15 

0.149 

X 

10 3 

15 

< 5 

X 

5 ) 

0.739 

X 

10 2 

15 

0.944 

X 

10 3 

15 

( 6 

X 

6 ) 

0.142 

X 

10 3 

15 

0.591 

X 

10 4 

14 

( 7 

X 

7 ) 

0.243 

X 

10 3 

15 

0.366 

X 

10 5 

i4 

( 8 

X 

8 ) 

0.384 

X 

IO 3 

15 - 

0.225 

X 

10 6 

13 

( 9 

X 

9 ) 

0.571 

X 

10 3 

15 

0.138 

X 

10 7 

13 

(10 

X 

10) 

0.812 

X 

10 3 

15 . 

0.840 

X 

10 7 

11 

(H 

X 

11) 

0.111 

X 

io h 

15 

0.510 

X 

10 8 

11 

i (12 

X 

12) 

0.148 

X 

io 4 

15 

0.309 

X 

io 9 . 

11 

(14 

X 

14) 

-=j* 

OJ 

• 

o 

X 

io 4 

15 

0.112 

X 

10 11 

9 

(18 

X 

18) 

0.545 

X 

io 4 

. 15 

0.145 

X 

10 14 

6 

(22 

X 

22) 

0.103 

X 

10 5 

15 

0.177 

X 

10 17 

2 ‘ v 

(26 

X 

26) 

0.179 

X 

10 5 

15 

0.907 

X 

10 18 

0** 


C 


IDGT is the approximate number of digits in the inverse which were 
left unchanged after iterative improvement. 

** There was no convergence in the iterative improvement. The upper 
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bound on the condition of Cj was 0.202 x 10 ^ based on the "best” 



Note : Subscripts R and Q on- IDGT refer to inversion of R and 

respectively. 



is computationally impossible. Based on condition numbers and accuracy, 
either formulation is acceptable for (N+2) - 5; while for (N+2) > 6, 
the first formulation is preferable. 

It should be noted that Finlayson (1972, p. 35) alludes to 
this problem, in his text but the comment is quite qualitative and some- 
what obscure: 

"The orthogonality of the polynomials gives 
computational advantages, although the same approxi- 
mation can be expressed in terms of powers of x, 
if the computations can be done accurately enough." 

The preceding analysis used the roots of the so-called shifted 
Legendre polynomials. These are defined by equation (45 ) if one lets 
a = 0, b - 1 and w(z) = 1. The polynomial coefficients were computed 
using the relationships of Villadsen (1970). Figure 6 shows the behavior 
of the first four of these polynomials. The roots were computed by 
shifting the abscissas from Gaussian quadrative formulae, available an 
Abramowitz and Segun (1965), Love (1966), and Stroud and Secrest (1966^. 

Although most of the problems solved by others using orthog- 
onal collocation have not required over 12 collocation points, the 
results of this section point out a computation disadvantage of the 
second fonnulation which appears at a fairly small degree of discreti- 
zation and which gets progressively worse. The first formulation 
requires some additional information but successfully circumvents the 
problems inherent in the second formulation. 





• C. Sample Application of Orthogonal Collocation 

The preceding can he applied to illustrate how a partial 
differential equation of the type encountered in chromatographic systems 
analysis can be solved. The same example will be used (except for the 
forcing, function) as was considered in the evaluation of the finite 
difference technique the problem is: 

(i/ Pe ) /Jz 2 ) - dv I di = 

y (z, o) ' = o 

y (o, o) = 0 (e) ; © > o 

b y (Z Q , «)/e) Z = 0 ; O > 0 

The preceding analysis has been conducted based on the spatial interval 
of as the interval of orthogonality for the orthogonal poly- 

nomials used in the trial function expansion. However, as was shown in 
Section IV, Z Q should be different from unity depending on the value 
of the Pecl~t number.. To avoid the derivation of addition polynomials 
orthogonal on an interval j^O, and the determination of the 

required roots, the above problem may be rescaled in the spatial domain 
by the following change in variable: 

Z ^ = (1/zJ z 

new w 0 

Therefore: 

Vd z = (1/Z 0 ) V^Znew (70) 

(i/jzf = (1/Z 0 ) 2 (i/^z new ) 2 
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Use of equations (TO) and ( 71 ) and deletion of the subscript "new" 

yields the re- scaled problem: 

2 

(p!) (gy) (<jV<)z 2 ) -^) 0y/<)z) = <^y/a« (T 2 ) 

y (z, o) = o (73) 


y (o, e) = jzi. (©) ; © > o 


dy (i> e)/dz = o ; e > o 


(7*0 

(75) 


Where one was concerned about the dimensionless length of unity in the 
old coordinate system, one is now concerned with the dimensionless 
length of (1 /Zq) which now corresponds to the outlet of the bed. 

One can now apply the matrices given in equations ( 58 ) and 
(, (59) to yield a set of coupled, ordinary differential equations. Since 

one requires the partial differential equation to be satisfied at the 
N collocation points, the result is N coupled ordinary differential 
equations. If 



i = 1, N+2 ; j=l, N+2 


( 76 ) 


then 



i=l 


0=2, N+l 


( 77 ) 



kQ. 


Application of the boundary conditions yield: 

y-L = y (z, ©) = f («) = 0 («) 

and 

N+2 

■ 0 = I O s _ 1 ] N+ 2 ,i. ^ 

i=l 

Using equation (79) one may solve for y N+2 (=y (Z^ +2 , -0)): 


N+l 


y 


N+2 


Zfa [ a f % 21 

[= = _N+2, N+2 C= =" Jn+2~ 


y 


l 


(79) 


(79) 


(8o) 


Equations ( 78 ) and (80) can further be used to reduce equation ( 77 ) to: 



; j = 2, N+l 


(81) 



Now, define the following quantities: 


Y = 


— — \ 


— 

y l 


y 2 

y 2 


y 3 

• 


■ 

• 


ft 

* 


. ft 

y N 


y N+l 

— 


— — 


(82) 


Y £ 


d Y-Vd e 


d . V d 6 


(83) 


Hat "D3a + i, i+ i - (XU, N+2 [a r% 2 , ij 

[S I ]j)+2, M+2 ; 1=1, N ; 3=1, N (8h) 

Hj = W J + l,l - H J+1 , N+2 [s S' 2 l + 2 ,l/ 


[s § _1 ]n + 2, 


N+2 ; j=l, N 


(85) 


Equations (82) through (85) may now replace equation (8l) 


by: 


Y = A Y + b 0 (e) 


( 86 ) 


Thus, one has reduced the distributed system to a lumped system via the 
spatial discretisation given by application of the orthogonal collocation 
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method. Equation (86) is a general form for a forced, linear system 
of ordinary differential equations . The stability of the system is 
determined solely by the characteristic values or eigenvalues of the 
system matrix A. 

An eigenanalysis -was made of A for Feclet numbers of 1, 

10, 100, 1000, and 10,000. The number of collocation points, N, was 
varied in the sequence 4, 8, 12, l6, 20, and 24. The shifted Legendre 
polynomials and roots were used in the analysis. The value of Zq was 
held at 2 for all cases. 

For the cases of Pe of 10, 100, 1000, and 10,000, the 
eigenanalysis yielded eigenvalues with negative real parts or the 
calculation was stable. For the cases of a Pe of 1 and all N , there 
was a least one eigenvalue with a positive real part, revealing an 
unstable computational method. 

While stability is indicated by the negative real parts of 
the eigenvalues, an oscillatory behavior was indicated by the presence 
of imaginary parts for a majority of the eigenvalues in each case. ' 

Prior to this eigenanalysis, some simulations had been made which exhibited 
damped oscillation in response to 0 (6) being a unit rectangular pulse. 
Additional simulations were performed after the eigenanalysis with 
similar oscillations noted in the responses at the collocation points. 
Figure 7 shows a simulation for W=4, Zq= 2 and Pe=10. The forcing 
function, 0 (-0), is a unit rectangular pulse of duration 5.0 dimen- 

sionless time units. The value of the response at Z = l/z^ is inter- 
polated by the formula: 

y (l/z 0 , O) = r T R -1 y 


( 87 ) 
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where : 



1 

(VZq) 

(l/z 0 ) (1 - 1/Z 0 ) P 0 (X/Z 0 ) 
(1/Z 0 ) (1 - 1/Z 0 ) (1/Z 0 ) 


(1/Z Q ) (1 - 1/Z 0 ) Pu.^l/Zg) 


(88) 


While the simulation appears adequate for long-time dynamics, 

it appears inadequate for short-time dynamic situations which are in the 

pulsed chromatographic system. 

It is concluded that the trial function proposed for work 

in axial-diffusion may he adequate for steady-state analysis and for 

long-time dynamic analysis of systems forced by inputs such as step 

functions. The next section presents ideas as to how a different trial 

C 

function expansion may give better results for the short-time dynamics 
prevalent in the pulsed chromatographic system. 

D. Use of Orthogonal Collocation in Chromatographic System 
Modelling 

The results of the previous section have forced considera- 
tion of a trial, function which may be more adequate for purposes of 
the dynamic analysis of the pulsed .chromatographic system. For the 
problem considered in section VI. C., an alternative trial function 
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expansion is proposed which may he more suitable for the pulsed system. 
It is: 

N 

y(z,o) = f(e-z) + g(e) z + z(i-z) ^ a ± (e) (z) 

i=l 

The difference from the previous trial function fom resides 
in the first term of the trial function, f (9-z). At Z=0, this is 
just f (-9) or 0 (■©), the input wave. At subsequent axial positions, 
the trial function is the translated input plus some additional terns 
to 11 correct" for the axial diffusion. This analysis corresponds to 
the situation where Pe = oo . A hyperbolic problem results with the 
analytic solution: 

y (z, q) = <t> (e - z) 

From the problem of diffusion and convection alone, this function 
included in the trial function should yield a response with a correct 
mean since the mean of the response is unaffected by the diffusive 
term. 

Use of this type of trial function revision is presently 
being studied. It is thought that results obtained by using this 
revised trial function will be more favorable than the results obtained 
using equation (46). 

VII. CONCLUSIONS AND FUTURE WORK 

This report has summarized work conducted during the period 
June 197 2 through August 1973* The work has dealt with the area of 
chromatographic systems study referred to as Model Improvement. 



Previous work dictated the formulation of a model which took into account 
more of the dynamically relevant transport mechanisms. A model has been 
formulated which includes intraparticle diffusion and rates of adsorp- 
tion that were heretofore neglected. The model has been analysed using 
the moment analysis technique. This analysis of the proposed Inter- 
Intraparticle Adsorption Model indicates that the gross characteristics 
of actual data are more adequately predicted than with previous models. 

The mathematical complexity of the proposed Inter- Intraparticle 
Adsorption Model has prompted consideration of numerical techniques 
appropriate for the solution of the partial differential equation models 
which are being postulated. The use of numerical techniques for the 
second-order models being considered requires the use of a finite 
terminal boundary condition. Criteria have been developed for a simple 
model' 'wherein a finite terminal boundary condition can be applied which 
yields system, responses which are for all intents and purposes equiva- 
lent to the responses obtained using an infinite column boundary condi- 
tion. • 

Investigations into the merits of the more popular Finite 
Difference Technique for solutions to partial differential equations 
of the type encountered in chromatographic modelling have concluded 
that their use is not warranted because of the large degree of spatial 
discretization required for numerical stability'. This drawback forces 
the use of a large amount of computer time to perform simulations. 
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Investigations into the merits of the method of Orthogonal 
Collocation as applied to the solution of partial differential equa- 
tions of the type encountered in chromatograph system models has 
resulted in some inadequate results . However, the results are based 
on a trial function expansion which is more suitable for steady- state 
system analysis. The use of a modified trial function should give 
results which better model the dynamics of the system; 

Future work in the area of chromatograph system modelling will 
deal with the establishment of a reliable and efficient technique for 
the numerical solution of the chromatograph model equations such that 
the dynamic effects of the added transport mechanisms, which prohibit 
analytic solutions, may be adequately studied and analysed. 
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■ NOMENCLATURE 


A s 

A 

A y ,A X; ,A G 
a 

a i(«) 

a ' - 

P 

b 

To 

C i ' ... " 

d.(0) 

d 


f(S) 


f 


g(0) 


unit impulse, Dirac delta function, 
system matrix . defined "by equation (84). 
areas under output response curve, input respons 
curve, and impulse response curve, respectively, 
lower bound of interval of orthogonality used in 
orthogonal polynomial definition, equation (4-5). 
time-varying coefficients in trial function 
expansion . 

ratio of interfacial area to packed volume, 
upper bound of interval of orthogonality used in 
orthogonal polynomial definition, equation (45). 
forcing function vector defined in equation (85) 
scale factor used in orthogonal polynomial 
definition, equation (45). 

time varying coefficients in trial function 

‘ . C 

expansion. 

vector of time-varying coefficients defined in 
equation (63). 

time -varying function in the trial function 
expansion. 

vector of time-vaxying coefficients defined in 
equation (50). 

time-varying function in the trial function 
expansion. 
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L 

M 

N 


N. 


HU 


N toG 

Pe 

Pe. 


Pe 


E 


P ± (Z) 


S 

St 


r 

R 

r a 

*1 

IU 


length of chromatograph column. 

equilibrium constant. 

number of collocation points. 

the number of reactor units, a dimensionless 

measure of the rate of adsorption. 

number of transfer units, dimensionless. 

Peclet number, dimensionless, 
intraparticle Peclet number, a dimensionless 
measure of diffusion rates within the particle, 
interparticle -Peclet number, a dimensionless 
measure of diffusion rates within the carrier gas. 
group of polynomials, initially arbitrary but 
later constrained to be orthogonal on interval 
[a, by equation ( 45 ) . 

matrix defined by equation (62). 
matrix defined by equation (66). 
matrix defined by equation (67). 
intraparticle space variable, dimensionless, 
vector defined by equation (88). 
particle radius, 
rate of sample adsorption. 

moles of fluid in particle per mole of adsorption 
sites. 

moles of fluid within the total bed per moles 
of adsorption sites within the bed. 
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R 

Rl 

R§ 

Res 

s 

W(Z) 


w 


y 



matrix defined by equation ( 49 ). 

matrix defined by equation (56). 

matrix defined by equation (57). 

residual formed by trial function substitution 

in a differential equation. 

Laplace transform variable. 

weighting function used in orthogonal polynomial 
defining equation ( 45 ). 

weighting function in weighted residual integral, 
equation ( 43 ). 

matrix defined by equation (76). 
adsorbed phase concentration, dimensionless, 
interparticle gas phase composition, dimensionless, 
intraparticle gas phase composition, dimensionless, 

equilibrium intraparticle gas phase composition, 
dimensionless . : 

vector defined by equation ( 48 ). 
vector defined by equation (82). 
vector defined by equation (83). 
axial position in column, dimensionless, 
collocation point or end point, dimensionless, 
axial position where finite terminal boundary- 
condition, equation (33)? is applied. 
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GREEK LETTERS 


P 

As 
A z 
8 ( ) 

€ 

■6 


n 



particle porosity or void fraction, 
time increment in finite difference method, 
space increment in finite difference method. 
Dirac delta function. 

Krone cher delta 

void fraction of the "bed. 

dimensionless time variable. 

the n [ moment about the origin defined by 

equation (l4). 

the ri^ 1 moment about defined by 

equation ( 17 ) . 

function which satisfied boundary condition 
in trial function expansion. 


SUBSCRIPTS 

I 

II 
i 


d 


refers to Case I boundary condition analysis.- 
refers to Case II boundary condition analysis, 
refers to space level in Finite Difference 
technique; refers to column in Orthogonal 
Collocation matrices. 

refers to time level in Finite Difference tech- 
niques; refers to row and/or collocation point 
in Orthogonal Collocation matrices. 


MISCELLANEOUS 


[ 1. 


03 - 


ref ers to the matrix element of the row and 
i th 


column. 
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Appendix A • Derivation of the Inter- Intraparticle Adsorption 
Model Equations. 
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/V 4 V &£ USEto: 

(dy*/dr) ~o ; -e>o ( xr* is) 

f-o J 

£. APPPbpp/A-rE -zrAfrze ppet/ccE 8 co/)Mev Ccfi&mbAt'&. 


Jim. uTi-t) = -F/nUe. ; t?° 


C*Xv) 

( txvl) 
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# APPZsPZiAre /Nival CoMbmotts 


y( Xj o) - o ' ; 

X 7 0 

CxyvM-) 

%Ch, r ) 0 )^o . 

. \tc> / r?o 

UcxS/m) 

t a CKj^p)^o ■ 

x ? r zo 

Ocx*r) 


-2 oEAjzZ to mtz- Pzeviafrc v &e£/ve/> amj> 

^057 VL-ATEfi £<$> O AT/ 0/05 MOfLE £,^AlE&A-L. ; XT' IS AJE^ES^/P£y 

/\/n> A?P£oPp/#7£ _ -TO }>£f=WE /trt/y iMpt~£/yij£A/T I>/M£A)SloNC£& 

VAE/A&L&S 4ajJ> ObEPE fC/£AlT%. 

OPE MA\{ £)££-/*}£ THE FO^covO/VAr QOAVVr/ES: 

A A J>r/* 7 e/V$fOAJt,£S<; t/lyiE : 

o = V-t/L (xw) 

l- 4 A>//neAl$/DA)i££s> AX/ Ac LEAJC TA ■ 

H ■=■ X /L- (XXYi) 

3 . 4 J> mBfefctoLESS ZAPial LEMS ?-# : - 

si - r/P 

4. A brty£A$fOpo£$s AVAL MFFOS/O/J M IM8&Z'- 

?e £ - V L-/i>£ (.XXXjlaU) 

5. A b //tIE/VS/O/VC ESS A14ZS TpAAXFEfZ CoEpE/0 /£AjT •' 

a f t- / a ^ { ir KX*v) 

£. 4 j>/ME/0&0Atce$s XAJTPAP4tr/CL£ J>/EPt£icjJ A turtlgEe' 

Pe A = \//-/A* (XKK^) 



7 7//r £*770 of rto ces of Fa/(D pH£sjev r jm rpB PAByidtE 
70 Tff£ Wolf's. oF A’bSoZP'r/ON S fT&s: LOfTA/0 T/+& PA £ /'C L B: 

R z = pp 3 /p L (i-p)C ar (xxxri) 

&■ a J>irt£fiJS/QA/C£S s &AT£ oF APSofcP TtO/V CcFFF/C/F/Un 

M ru , 85 k«<k L /v (kxt>b*) 

9 A} 2>fA7£/V$f0AIL£SS AASOZBED PaAsE do# d£ASTP/f-776A / V 

Xq, “ C*. (■ x k x ) 


T//£ P/ZFCFP//V6 0£P/AJ /r/6WS WAY BE APPC 70 ru £ 

bffMeius/OiU Al sf r of £<?(/4T70a/s &fUEAJ. 4FT£fe tF£ 
5o&$£C}0£A)T ,7?Aaj /PUC4-T/CAJS , A S>&7~ OF 2> frt7 e/JSi'Q/U LSs»^ 
FqvAT/otQs 75 Z. &Af/Z£j>. T//JS PE&ULTS ace: 


TF£ -Z/V Te£ pA&rtcuz a. As P44^E M 4 ss BAlAajce A : 

^ jdl .A,, ru-n-) = dy (KXK^x) 

<J& 


t£ 


dz 2 - dz 


**< (y^Lj 


7 F£ FajTZA PAen<U£ ■ CAs A 0 ASE Fl^sS? i* 3 AFAFC& /S : 

t (iTL Si -pf) -w/j* w 


the AAscfzesp ?fa$£ mas-s %Ase .* 

c*y) 

A 7 FF£mo£>t//<sa/? 7 /C FBLA-TfOA) SASp BETU/££A Thf£ =»TRAPHericL£ 
Affb 40Sogp£d PaAs££ : 
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M £xre&v#L Cx.fi/rez- x^a/tza P4£t(ccb) gousud/HZV Coit&moA/ : 

= A/£o & (y- y^); 

- jl ... ( yJ<U+) 

4-V xp reZAi/lc C *pt/ 24 MzUccx) Boo/uoaeu co»z> fT /ou : 

dy*/d/z, = o j /v ~ o (jtJ*\/) 

4PP£cP£(4r£ XP7£Z PA&nctg g<DOA)t>AP{ c*4Jj>/rz oajS! 

=■ 

M*k, » <9>o (xgi/*) 

£ — ? 

ApP2oPZt4T£ XAJfr/AL &>//£> trsoAJS' 

yc*,o) =o • dtfvju) 

% CtyKjO) - a ) Cy^i/^c) 

Xa($ y sZ y o) - o ; Z-?o,/lZo .(yjfiCX.) 

/7 ~ 5#oOCj> B£L A/OTXb 7t/4r PPXSE.KfTA-'TJOA) o R ABqVIX 
equfrrt oa/s //v rhfe rta/A/ £o&y of TU£ £xpo£_r may om/r 
JH£ S<J&SC£/PT o. op X-o. A/[jJ> 'THE P4&EPTU£sf SE/) 

fofrsce/fir te) oaj AS C6 ) * 
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fi)DfliefiJ&4Tv££ 

<2p - £A7 /c> 07 = &HreeFAc/4c 4£eA to wc*:^ volume, C' 

Ajft)= t/A/fT /rtA/cs E F*>Z R&AC -7//HS. J 2>/ZAc- PectA TlHOT/oK/, 

44#) - uurr /a/Puc^e fo/l P/me/vs/oFcEss 7 /ME , d/Enc pelta FoKxETibH. 
c A » motes &f 4z>soe&zj> gas pec u/u/r wettur of a&so/z&fa/t. 

C*T * MolES OF HbsoZ&Eb £4S PEE OAJ/T U/E/LHT OF AbSoQBEkT 

- «? HE/J 4 ll 4&L>ofLf>T/oH s rres AQE occuP/£i)- 
J> A - XAJTfcA PAfCTtCLE D/FFVS/oAj COEFFICIENT 
2>s * ^AJTeZPAVThllE. 4X/FE d/FFUS/OAJ Co £ FF HLi£A}T . 

Gr =; OlocAZ vecoc/rv of gas z>T£-BA/y). 

Crg — MOC4& VECoC/7 - 1 y OF* c A-fLF/EfL &/4S. 

It *4t g =• EbLofLpT/ \j£ g-ATZ Co/VSrrtfi/T bEFWEfi as TUe T//HE ft ATE OF 

c hA/jce of the of /Holes oF AbsoZBsd fha^e Ae£ 

V/JlT VOLUME* of PAE.7~1C.CE TO THE /HOLES' OF jr/UTEA PAPT/CtE 
PHASE/ PE£ OA//T VOLUME OF PAFTCLE Vo /A VOL UME. 

ArVj - F/c/r) T&WSFE/Z. CoEFF/C/EAjr bEFWEb AS the F/lM E4TE 

OF TZAA/SFEtc OF Moles Pe<£ EKTEfCAjAE su&f/KE AP-esA 

Pee pazt/cee, 

U - LEjv<S7H of eFeon/A-TtxZEApA C-OLUm/V 

/n - eyo/L/Se/v/v Go/js taa/T AeF/A/eA /aj ecpuAT'oA/ ( 1-0**)- 

/v 1 ~ ecpu/L/SCfoM CoAJSrAfijT b£F/NEb /Aj EQVAt/oAj (juU)* 

=■ A/VM6EZ of csPAEPES P£re (JAflT VOcV/YlE OF & EJ> . 
f/ffa = bim£/VS/0/VL£s,S A0S0PPT76/U (LATE CoEFF /C / EA/77 
~ b/MEA/S/O JULEPS /HASS TPA/VSFEfl CoEFF/c /E/vT* 
f p d PPEssuZE TAJ TAE &£D. 

% - J>)/hFAJ$/0H CESS 3z/JTKAf>A£ncc£ b/FFu^/oA Coeff/C/e/J F 

Pee. ~ & f meF £ t o A less jCA/ TeP PAqtic l b Ax /Ac. H/ffvs/oA] <Zoeffk.i£A)T 

T - PlMEFS/oA/pL £aO/AC D/STMA'CE WITH/AJ /he PAET/C-UE. 

r a* P/M FA 'S/Ofi/A-C POS/Ttb/V UHTF/fj THE PAQ-Tlc.CE ElE/VEA/T Op 

th/c/c pess a r. 

71, - P/M EAJS/O/VCESS -T/JT&AP AfFn C-LE LEfJFTU- 

R. = PAfLTlCUE FAA/VS . 

£$ ~ GAS CctOs tAaj t. 

Ej: •= e AHO oF Moles of FLO/P oormw /HE PAF/TCCF to THE: 

ftiOCES of 4b30£pr/0AJ 3>/TES /A/ /HE P4£ Tl<LL£S . 

$ = CE OSS -SECT/O// 4 L Blew EFT of tFE BeA. 

~t - P/MH/VS/0F4L T/rtE . 

~T = TbA/PeFA tu£E of tee 8 fa. 

V = L(AJ£4E /AJ7eesr777AL l/ECOC/TY. 

X " £>/E}£A/S/OPAc flt/AL- PO3/T/0/V /# '/EEL 8E&. 

y_ - filMBNStoUAE A//4L Po$/T/oA/ to/TH/fj THE Ft£fY\£/[}T A/ . 
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Xc* = r*ioce: rZAcrtati or picked 4^so^'noKJ suites, -dumenSiouless. 

<£ =• coMce^Te^vr\o^ im the ^u*reej>A.e^<u-E. ^cv&w. ) fciittEKisao wL ^ ss ' 

^ . " C^HCENTtA-TlCi^ \k* THE 3:NTVZJVP^ei-vccE. J2.£G\tifJ , JHtt\£V UNLESS.' 
vjj = EQO U- vB>£a(J tn eCs^CttiT^^n 0*J iu WE AVJT^A PA k-YV Ct£ 
fc>\IY\E*aS\&fc}L£ss. 

2 - b\rvvewsvC^L£SS t^CKjr£2^AfcT'<LJ»E-') LaNGm 

&|2£.£k: UCTtees 

£ - pfV^r^cL-E- Pofto^aM op El2AcYt<svJ. 

6X ~ Pvnr\EfOsmKiAU A-f { A-L. Yjoo^g^&siT umanu T-M-e fe&D. 

<4 r ~ £> tmjaos ioual p-adwi-l- XMcuz^n emt- with-uo the. Pp^xxqxj^. 

6 « ^ELb Vovt) T^ACTTCjW. 

= tvioLAR i>E^lTH CsF THE- G»AS. 

■fu ^ APPAK-E^t ;i>£KiS lYM OPTHE • Pfye.TVc.UE y>j£^H-T oP 

Absoe-&£NjT- o^vr voLomE. of paethcile* 

© n MraeNSimLE-SS TTflse.. 
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Appendix B Details of Moment Analysis for Pulsed Systems. 

details of momfajt ahachsie fq£ polled s ys te/ps^ 


~T(te IMPULSE- FESPo/CSB OF a OHEM>t*F zeaotcR- /% SCSHFT/MFS 
Z.£-FE£Z B-& 7 0 A* TUB CES /PEACE TfrtlE FFEq U£Al Cy DlSTtlSUr/oA} . 
THE STEP FEsPoHSE : } THE iHT£Q.pAE OF THE /MpUL S £. /Z&SpoH S B > 
XS Tff-EA T£-Fft\£-0 THE CUM UL-A-Tl V£ R.E S ib£/J C£ T/AlE D/S TJL t£(J T/oU* 
' THESE $u4-A/rrr/£S ■ EBSEMgLE THE J> /^r£/Su 770 AlS STt/D/FP //tl 

statsst/cA-c AH/feys/s • x+fi. f the PEogAS/oiry o/srEttSurioAi fohct/oh 
A tob 7 7+E & o /n V L AT l !/£ b/sTJZ/3 i/r/OH Ft) Aid T/alf. 

ONE Af AY EAAfDACTjZfZlZE- 4 M&TICULA& D/sTfLf BUT/CA/ 

VOtJCTtoF) -Fee); BY JITS MomzA/TS. , tAe MOAIEKSTS CA/U &E DE-PI A/zd 
BY THE Fo^ow/H^. £<3uAVotJ; 

«v? 


THE AfCWEA/rs DEPWEA SY THE H&oYE &qUA-T7oA] /4EE 7 E £*?££> 
4$$0 Cut£ A7O07/ZATJ OF A10 AlBA TS A&OVT THE oElG/Aj, THE 

b&uonu/Vjprc/z. /a / e^u&noA/ 6 %) /$, Tugr the AFba c/tuhEZ the 
C oFi/£ J Fco) . 


THE Ee LA7ioAl^M/PS OF THE uAZ/ous Aiom&A/te To THE 

LAPlAoE TfrAAispoCA? C-AH EE J)/S HELMED \MHEtU T^-P tEAH^FdFA] 
is iHE/npH/ /aj te&ais op A/u T/jp/v/rE series. 


LET : 

THBFE FO££! 


■fCs) = ^ [_~Fc&)] 


«o 

-Fc$) == s * -fee) do — J ( /- s& +£ 

O O 

CAeev/HC 7H£ /A/retFAt thFho&h aa/a f-ActoF/IQL out J-f&de ■ 
6 tups: ^ ° 


?(s) = J fade U- 


sje-fcdb P*fy Je 

J? + Fs w . 


fFfe) do 




-t-s* 
- 6 s . 


fe* 4 te)d 0 


j_ 4 fe^lbk/s 

+ 11° 


7 -Hid* 


He so /he 4 s- peers op th/s 2 >£uzLopfo£fisT Pacha le^ those PEzsemtEL 
SV DOUg, las j T.M ./ 1 Process by/iaM'cs a'ld Q>n'fo/ } VoU*iet. : AA^t^sis oF 
Dyna/tuc System ” Pr<zA& eg,- Ha//; T/rc., Engfeivooct &<'fk;AJeu' T ^cy ; 
*472, pp ■ HSC-m6Hj 26S-ZCF. 



1V1 ■ j - fl. 
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ZgCAU-tNt Tfr£ DBFtMrrroA) OF Tfr£. Atf* AtO#1£Ajr C/f/g/J /AJ 
EtyJfiTtOH C*-) j T/f£ d+fij MOW &£ Ltf£tTT£A/ t A t TEgfrlS 

OF T*£ /Ytortie-A/Tl : 

= (jrft f 0 ) Js ) l 1 ~ + i s’-zu^ - 1 + J- 4 -s*ai 4 - • ■ -J 

c)/</^ CAt v sfroia Tfrs. &££AT(o Ns.fr ip of Tfre: c.AP^c£ 

rBAA/SFOE/y} TO t 4£- frfo/H£JUr FUAJC,r/o^ 

«o *? 

JUm J — J f to = 4 -f 

S-» Cc ° o 

*»#£££.• A f = TfrE A ££fr UA/&££ Tfr£ Cv£V£ 


4 -M*' 


Jam t- f - JUm 1 &<z SB -f&)de> 

<2> 

. •'• h C. = -AfM, 


ALSO 


<K> 


jNm €. [ fe^-foJs = 
s+o as z K J 




jUr, Jo Z £~ S * ffe)d& 

5-36 O 

CO 

’= /© *f(e)Je 

o 

= AfMz. 


-TfrEse £ Xteess/otic. ■ u£/& -TO TfrE Mc£.£ C£AJ££A£ Z£CNV<Jtefr,fr , 

M 

£ r K £ ®f^)de] - c-iy jv-fea® 

5*>o ^ _i *> 


oe*. 


iwt — 0 I f 0)1 =■ L-if 

5-5o 6S° L J 



5o<-t/<wG FOE. /U n AND iMS£CT!W6 £>BF tfO l^LOfJ & F A-f 

YfELLDS *Tftt C_H liVt\c»nsyhp Foe. the. mowent^ twc 

t*VPcrtC£ Tfc^SFOt^’. 

/V 

( AC.) 

tf-?0 


M n - c- 1"? 1_ 

r» ->o dks 0 


L**V 


Tf»E. Pp£d£I>\N<Sr- Bt {/;&£/> 70 Ofgr/f//V 2V&>e/7 l 7/jT7o V FO/2. 

SV&T&MS Fo£.C&£> SV ACB/7Z4?C y pvC^£i /aPotS.. CoAfcs&£:£. T tf£ 
FoCC-CUJ/AJ^ -TZstHsrZfZ. FVUCTtOfi/ &C.OCMC A/ti&g* 107: 


Vcs) ^ 


WVEp£! 

Lxte^ i XteO IST-«E PofcCt ^ G. FOrten&M, 

L^Ce')~] \ ^(- & ) is "ruE. output vaj^\i£. 

SH$T£«\ TeftHSFEt FvJMCTlOK), UjV\CW VS 

t«e Laplace tplavj5^q£t^ op t ne v^vol^E- 
fc€s*to^sE op nue s'\ , ste.«c\. 

TrtL g>toc_t< DlA&Elft'm y l£LDS THE Ftt-LOUAlftfr 5LeuarrioV^tt\P • 

Nl/O - y-Csb C^c) '' 

Nou/ t4<^A 'rEAV^PoEfia OAK) BE CLotv^ibCEE-b \Kj ~r^P_mS OF 

rr& P esP E crt-\ vje. m oWEM*ts - 


Yc*.^ - 
= 

£cs} - 


xcs> 


G-Cs) 


*<*>* A ? [1- 5 A,* -^X* +U* 4 *** ^ 

YES') “ A^Li * SA. V isV^ 3y -F^4-^/V4 y-...] 

= Aq. [^ X - S/U,ft. + is l /Jz 4 - "fo -*- • * J Cv/a) 

SOBSatTOTlWG £ qoA-Tio^S (>i>0 /V) AA;b GfX ) XWTo EG?OA T toW i 


YIELDS TW£ fOU-OOJtKJfi : 
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><«, L » - -i s J /M iv y - ...] 

= -McU ' s C/tf.x+^s.) -+ i^CA^A/U^Mfc, +/U Z(r ) 

-i **^»x + +$Mz^, )l ’‘>7j Cr ) 

+ i+S«0u 4)< +4 /V, 6 M ilc +t‘// z&/ a Z j' + ^/U 3& ^, x 

- • • ■ 3 ... cvjU) 

E-<}04tvU& the CoEFficiEMTS OF THE POVUEP? OF -s> IN 
ecj^fr'noKj {Vjil) N »Ei-t>s Trt-e PqLUsW iM& e£LA-r<c>NSHt-K - 

4>j - A x A(, ^ 

4//y — Mty (six) 

Mz y " ^2y * "/• // 2< i ( X) 

- ^a* ^ 3^/c^zjc * 

E^tffVUofj CV/U*') ST^TfS ttt4T 7W£ PfLOhOcTT OF THE A£E/f V A/AE&- 
-tUE. IMPULSE ^E&P oUSE At,b THE AKEA vA/£>eZ- THE /AJpUr pUoSE 

Equals the AeE?A vmae/z the ootpot Fvc&e. 

equation (ax) st Ates that the HiE/fU/ of the output oc.cu£s at 
TUE isu/yj OF THE A7£AMS OF ThE //UpOT A NO) SL/y?Puc.SE. 


TH/$ DISCUSS 10 AJ Cou LA &£ CO/UTSA/UEJ) / A/ 7 EBAJE oF THE' A/N3&E& 
4&SOEU7E nioyyf£A)TS &ur fte n ?2 } jz/ute/zest ceajteps *4/ 
MOM £ JUTS ABOUT the TYJEM// j M, . THESE Aion?£A/T$ /HLE 2>ELp/VEJ) 

By -THE Follow no 6 • 


<o 


F„ - S (e ' 







n Z2. 


<X*L) 


EyPAHS/cH OF tHE A&6VF y/£U >* S EXPfZE SS/o/OS To fS 4C n JAJ 

TE&flJS oF THE MOJTJBA/XS M n • 

FoE a =. sl , tHE EEsoc.t is: 

Mz. =- 


Ff 2 . -* A( f 





M-l T4E \SAeiA4[te TAP PeynztfuTloM PUVCTioN. 

Fo£ n -* t 7 -fife zepvlt /s •• 

a $ ® At % - 3/«/^z +ZM,* 

VSfA/c> Et?v47.fo /V fXA*o) TAf/S jZEbVCES 70' 

M j * ^ - Ft* ( Xkv) 

i/s/A/6 eqvA-noA/s (x<U) 4 a/£ Cx*V) } o/ve cap eelueite 
EQ uA-T/OA)S (X) 4AJ6 CX4A) //O tEPWZ oP THE BE^PECTWE Md/VEa/T* 


AZoUT THE MFAfV j *t 3y , 
the peso fts 4/ca: 


; >^2 y , ; PMfysfvC EQOATfOrt UtJj 

Af z Y 

= Fx 

*■ ^6. 

c w) 

AFP 




M 5y 

r *3* 


c V?a) 


THUS THE ISA&/4A/CE. OF THE OUTPUT /■$ THE SO/77 OF THE. 
VA&I4MCE. OF THE /A/i>uT A^HF THE //*1/>C>P$E E/zspoAJSE . THE 
Tp/ZJ> OloryjE/UT AQooT THE OZ/C/AJ op tHE EE^Poa/^jE C&£C4-re£ 

TO SkzAU/) A? THE son 7 OF THE 7A/^P H10H1PA/7S A A OUT THE EESPe<=T/iE 
PI a AH $ OF' 7HE /SPOT p/i/6 XHlP^SE ££SPqA)SE. 

e 



Appendix C Derivation of the Transfer Function for the. 
. Inter- Intrap article Adsorption Model. 
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O£WV 4 T/ 0 // OF Fc/A)cr/0/V Foe Tff£ X HT£ £ -XA/T£4 - 

Particle 4j>jso£t>Tt&/>/ moo&i. 

th-e Mob&c eqUA-Tious a&£ : 


J_ __ C$9 

Pe e d2 L 


'^boc " %a)/l = t 


~ tl 



+ £%] c y «-%«*) = 

°>y+ 

d& 

C* it*) 



C 4sU- ) 

2* = ^ ** 


(*>/) 

[(?/*)**(£)//>«] ^ 

jz,~ d 


zz O > SL z O 

J 


(l/A.) 

y ( %,(>) ss 0 ; 0 


(VjU) 

^ A ^ 

X« = <? J, 2 


(v/tiL) 

C*k) 

y <: tf, * As - >4* 

^^71 &) ~ firlrf* 


(*) 

<*ji) 


Ol/E rvtA Y AjOoU PfaOCBEJ* OUST# T ^ £ f ot/CTtoAl 

OFfUV/ms/J: 



8o. 


&t£&Tj T/f-jcZ Ttf£ LAPcfiCB TMIUSPOZM o£ £q>UA7TOA/ fe); 
TV/S y/£LPS : 


■=» 3y fys) - 2C*>o) 

l/£w& £G?u4noA/ (\oU>) } TV/s Becomes.: 

t ^’ -^ • V»l„ - ^S'V) 


(X*Q 


you/ } r/ftc/z v#e cAP^-Ace T&AA/sea/zm oa zqvatioa/ (a*-) • 

t ( *T[ 'i *^]-M»[)i‘iv>-V?ow>] 

V&/yG ewnotUvtii) ) tv/s Becomes: 

4 (vYL d ^ T^* k 


ouu.) 


7£4vsfoe.m/A}& CQu/f-TtotJ (a*C) y/£cJ*s ' 

-KC W )] - A/eu Lyx(w,s)' %?(*,*, s)] 

US/A/6 £C?U4-770AJ Uk), ~ r H(S g£CCMe £•' 

£ & x* £$,€>*) - AfaCy+C%*,s) - 3*Q4s)] new) 

TT/JA/S PORJ71W6, EQV/pr/ofc (aV) V/£t AS / 

y/.dXjs) = />7 ^ ?,/?, s oV 

USWi EQO/frroto C/w) fA/ EQOAT/oAJ (XaV) C /(/££- 

~L^ s ytat^s) = Me*.[y/C%4s) ~ qtcs.'l.s) 

Solo/A/C* FgP~ $£(■%%*) G /l /£&.' 
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Now eqv&voA/ Cxv*) fh4i BE $v8swvTEb j:attd £<?t/&notJ LkaaL) ~ tq 
£ LimiuA 7E yfcZj/i/s ) . BFreJL Sorne. Ze A EBjW C£ m£MTj r«& So 8S7/ TO - 
Ttot/ ice.soc.rs /a> the Follow i N 6 EQU/vrioKj ; 

t'fj.?**! + 2=. d y^C-Z/L^s ) __ # u^c&*jS) - o CxV* U) 


WHSP-E: 

<z t — -fu^ot fa? o"F s 





vow e.Q{//Pnov ft #46) ntt/$T &e solved foe ^C^a.s). zeu/B/re 
EQOtlTtoM QctfJU) A&f' V 


si, d *y* 


+ J Jfcte.'cs) 

J/L 




CXsCx) 


Propose /f jolvt/oaj of rue Foznu 

C*,A,s) = tosh CfcC Z) + Sjz smh 

SI. /t» 

TUE Co*/sT*/tsrs CwtTU HGc.pE.cr To A-) MV%T BE Z)£7££W/A/£A 

4fPc/c./tnoA/ e>F me rZAA&eozme-b &ooAJd/tzy CoA/btlro^ 

£(? vA-tto*/ CM ) c /ires ' 

4 5=5 0 

4PPctc4rroAJ of 'the Boundary totub/rro/i/ E^>u4Trof\)co). 

fiT/t/^S (AF-ree /1)/tAStPVi./l'77t\): 

C* « E y£ s *,s)/[cb~Os/nh(ti%) + ^ CosJr C 0 ,)J 
loueze-’ 

b = tit* c j ( L fc) ^IPe^ (xx) 

T/feZEXoZEf 

y ; = y (i ' A ' s )j[(6-j)s,*/i(rs,) twA^u^i f . | ^ Cn ^ 

e?u4T/o/0()(X*) /S WE fiJEC-EESAEV CoOpC/AJG ZECAT/o/JsAlp AJ&E&E 0 

For solut/o/Z c>f EsQuA-TtoN CxiL ) . ZECAll EQvA-y/oAJ CXlL): 

L ^ ~ **“■ &*' s) " #*'*'4--/ e s 2 C *> S) lM 



82 . 


E</4to/$nAJ(> SQvArro/v CxxJ.) /ft /L-l Vi£t6S / 


%u,j,s) = <jC%s){ L 

A/OWy 


l Cb-/)s(f}A}(^) + \Ia, 


) 


t£T 


Its) ^ 


h sthhl0.) 


La-0 tinhtfa,) + \J5) Cos/j/ffij] 


( kxJU) 


'7V£P£ rro&£: 


y+c^z) = la) 


SV6£7'1TV7'/0AJ OF Ttf/S /A/ To £ GOAT7G// lX*C) #/J d 

£ l\t£S : 



Z&A tztADc* emeu r 

( XXAAi-) 


u/ftetLe: 

W = lt*>l + *} (xxW 

Xk; t#£ 3? j BQHifTfOA/ (XyAM, ) is A C-tA/BAfg , f/orwoeeA/eouSj 

SB Co A/b-oe.6 BZj Co /US TKAj T CoBFF/C /EA/T j c£J> tA/A-g V E> S^fSL/leSJ T/4L- 

SQo/bTtoW. -tab <?#A£Ac,r£/e/sT/c_ eql/a-tioaj os £-<? vA-r-/o/V C x/Ut ) 

(itj p * ' p - ~ ° 

T/t£. g-ooTS AELS! 

t 

T> = ± ±Mt 
a £ - - 


Ttf£ jocvt/oaj ~r& EyuAnoA/ CKVAXl) as of tab /=oejn: 

%(?,$) - <V etplVe] +c z *vp[P z *J CXxv) 

UjMBZE 

P, = Pe g O+ V l *- (4 LCO/Pcsi ) )/z 
/ z = ^ <f/- Y^-^rcs)^) ) h~ 


A/oh/ 17 iS A/E CESS A-fZ / TO DETEgMlA/E Tf/E COEFF/c KZA/7S. , C f 

AA/£> <1^ 



(xxvp 

C YxV/ti.) 

A&P&Y/Art equ* 770 tf (WM-) TO Eq>u*T*o/V C/tf) 'f/EtfiS f 

c, ~a 

#PPC V tA)& £c?VAi-(otj QX Va) TO yff-e ZB-bVCG-b £q>t/Ant>HO.u) Y/HZJiJ 

7f/^s T//£ rPAMZf&esn Fo& 7%T UA/CT f/nPOcSB e^Poajse G/Z 

T£AV$F£& POA/cnc>Aj p^. T7HZ. /A/T&e.PA&r-'C-LE; S&MAI/SO/t. /9-A-S 

Foo/ubj L.e.j 

yz>s) = trpLPz* ] 

jtt fs sorr*e>L£ ~rc> ujotz/c uj/th ta? tZaAszfb/z avvctioN f=o(z 
TH& CHZArt7ATOG>e.AP4 OUTC.B 7" AT l/Tii-f Zr&TlOAJ OF TFIS 

ytscDS sumnAfry) : 

yu.s) = e'fpZP^iQ 

ij, JI C £> XT * 

** ~ Sri ^(P^b 14) + feVpe e " 

fos> - [to-tot (/-}CS)) +s] 

b s/ohClci.) fL^b-A s/nhOR*) + 1 &, CosA (ift,)J 
A - »*><./[&/&) ( L te) «r L /&*] 


EQVA-TtOA/ (t) 7 £A*j$?=ofrm$> to: 

tyCo,*') - 1 j (rAfriA/C A = l) 
£<?*/ At/ o a/ Cx + i ) t fr/j/us Foz/hS t s » ; 
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Appendix D Derivation of the Moments of the impulse Response of 
the Inter- Intraparticle Adsorption Model. 

J>£ZJV4rioH fyZ TfrE MoMZkJTS dp THE IMPULSE. sg n /r 

T#£ - scA/re/fPP-zn<LCE ^soUPnoti tft&£L 


T#B L/tPlticE ye/tAteBoK-fn zocO jit>A/ FofE THE /*jPUt~*£ 
pEsPofJse o£~ t£-WSf££ FUQ£T7caJ of THZ /hotel, : 


LjO,s) - Zip(K) 

lOfiBZEt 

p x s (fk e j7^) - 

*<*> = [*&CrO-A<A) +S] 

Ms) ~ & s /nh(yq,)J [ (£~O s wh(VQ / ) + Cosh(\Jafy 


R,(s) ~ J 

a 

-Mtu. + y|/e„ + 5 

'(ip**? 

f 

0 s + Alfyfrjfiz) 

— *- 




Cti*) 
Otic . ) 
<m) 
(v) 
Cv/i) 


A- C ALCUL/pnotJ OP THE /VE/W , //, . 



. CvM ) 


jJn LtitljS)] - -£ 
o 


/t '°‘ v /r /£ /oec.ess^v to p£7#£mwz the /vvm&p&roe 

up z^vPr/o/y (yfcC). ^S/ajc €.q>of\Tit>H &,)■’ 


dltiCi,s)J = ewfa)'d fa) '(&) 

£/•> ° C/S. 

jLA>UJ &/*/& 7tftL OH AIL /> ~ £>£P&AJ0£AsT‘ TUAJC-TlortS 

lEPt/JE-j) 4-$ P/f£7 of- ijS) : 



85 . 


$! X 9 &P£GTtt>fi/ : 

J*n. <?, (*> - o 

S -9 0 

1 / 4 / A/ 6 . £'MOf> / 7 /}L's £OL& : 

Jhtt ^is") = i. 

£y J/vsPect/oa/ 

£0 * O 

£-»0 

4/0 

^ £ • * O 

if-** c 


a) 


Ou j 


exit') 

cyxU'J 


$i VZ//V 6 • 1 he Zet^noM T/f/fr E//r?jr oe 4 - PEopv<lt /s 

THE P/Zol>Uc T OE T/f£ C.//H/TS; £<?u4T70M ( Va*') C/fA/ ££. ZODUC.EO 
To: 


0 &//*£ Et$ 0 fiT( 0 fj (}MS)- C/tfES 

* = ..' -it ^ (x ^ 

A/OUZ 






JLr>i d_ £(s) 

X*o <XS 


atv) 


Mi* <lfa 

^ o ^ 


^^7 c/ /Z e $ f y_ _£ 

*Te> cTS 


(XVxj 


€(/#/// £uc£ J>/&EE.&EAJTtA-ncAJ 0 P E<^J4TioA/ £a,\/) £/(/£$ /f Oor>if>C/C/frEJ> 

erPEessrort , t#£- EsM/rx e>f= the <U>/n&//r*f?roA& 0/ cr uypEAPoc.Pc_ 

f& AJC-7T /f-EE EclJfljQ £ v E tPE/Ufo /fit & T7fE tf ! // > ££-&OC- /<^„ EO/*C/Ta/t/f 
«/ poU/EZ S&e/ES. APO JJp&P T4&A/6 T*f& <~/P//T. AfETEp. Sty 

M 4 V(f>ttcA 7 iotU the E&SOC.T /&: 
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t = - ~g (. t-t '/me jr )^/S Pe 't (x uU-) 

040c. SoBsrn otioU or rh/s jtmto Ttf£ PJ?£C£6/A)6- EtpupftOAiS 
h£U>$ (AFTE&. £ooi£ /vaajiPv cation) •• 


M / * * + 


j e^do 6/t/g-//f6 Tty/fr 


0- 6)fl (Jf i/aiZx) 

6 


0-6) 0 JL 
— 2 £* 


= K, 


o j OWES TH£ F/AJAC ££SCJ<LT . 


M t — (h&)P/6 +: *~/ry)& X 


(Kiaaac ) 


6. C/tc-Ci/csfr/nAJ or the *$aao/v2> mome/jr 4 #qut tub 


nt£Av , ac 7 . 


l/Sf*Jt> £ <?vA7,cA/ T#e £**£&$ to 0 F'OK a 4^ IS t 

C 


Mz. ~ J^ru ^ L M/ 

SHPO ^ 


5-9 0 ^ 

i> /ffeiev t/a n o /u y/f CJfS t 




^ 

pgeroan/vc, tv? /A) a/cate: g> pyrAejz&Af rrAVcA aaj& oc/a/c tag. 
UA7tr$ «Sft/CH Al*££APV Be£A/ FooMP cnsgSi 


//& ~~ 2-/4/ 

*c 




s-*o ds*~ 

?e pFae-At /AJl 7 He ///£> icA-reo bfFFBpertT/Ar/oN QfVBS: 
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P e ^ s^o ds* 

P£PPoZtN(, TVE PlfFEtzuriKT/od **> T*Em T#E L<M/rS> Of TVE 
Z&^UCT /AJVOLV tU(> HVPEEBbUC EUAJCf/OAjS [&S(A)& POWEP &Z£/£S 
&!t /£$ ; 



VOiV 


dli 

as 


_ f M £u Z '*r &r- 


Cs+ ^f,/> 7 /^-) 



4 /ui> 

J**r _ I-ZAttltotr 1/&Y& 

as* L <*+ *&*,&)* J ^ 4 


Ock) 


(ju^) 


TAE//VC the c/sntrs Equ/frtOAys C*k) Avt (tit), svEsr trartotij 

zuro zqu#nw Oobt )j Hf't) some /vAA/zEv^AT/e/v cive<s. tme g^suETz 


*l£ * .ff(s.)% * N e '0 \ + 

^ 6 t W ly-' /S 


Cf < A 6 .Cul 47 tbN OF THtr THlRt, Mome-AtT ASoUT T^f OJEA,d , M? 


Vi, (A 6 EyumON tVA*i) Th 4r EtfeesStoN fog- /?) /s : 

/S } - -JE <EL " 3 M, M z - /if 

s ?° ds* 

?££foert/A/t> THE c> / PfEg E/vTtAT/ol) Afljj) s/ fit At / EYjajC, somEU/EAT 
6 tl/£Si 



f£&71 t E/S £ yfBEfS/oAJ) JH-E PeotEEUZE /S / C A f*&V4lQ . J Bur 

TEP/OUS- EqUA-TtOAj (*X**l) ' f/EsT £e.AJcE.£> To; 

yt/j = C /U /Et z Jlf+n, d* % *S j 

f -9 o f/s ^ 



IX#**) 
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Tft£/J To : 


// 2 “ 


'^0 4 - 4 *, c f^J p _f. Q>M / yU z 


P&£to&A)/A}& TUB l>iB/=£ e&Ort/lT/bAf TAkrtNd THB Limits Of TH£ 

Zbzult /mtLi /!*& tfvfieesoi/t* toa/ctioms tvs/vc P*w££ s£&£s ^PwsmMS) 

&W£$: 


^ - 


+ b&G J**tA /* +2i) 

Pc £ 36 Sj>0 &S 3 L> ■5**o <^s l/-r *3 toj 

frftXr fjb+n <^4/ ) / V<? di.-l- — .2- — ^ /••I 

" \5->d <#y ( jrs* ^5^ ? 4*/ Oo^O 




{x*t/) 


TAt/A/c thc c<Aitr &p e$vp-nc>v Our) > 5i/0$7/ rvrwi tr jM7t> F/pi/pr/on/ 

flrXW) ) Pd£l/tov$ fL&zoUTS) AyvP SH4/J iPULAr/Ctf )T/eU>s T/JB Peso IT: 
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Appendix E .Case I and Case II Transfer Function Derivations 
■ and Moment Derivations. 

C./}S£ X AA/J> CAse U TZ-AAlSFEZ FUNCTION J>EClVH-TloMS 
MO fYIENT X> t.£ tV^T /oXsT 



A- 


CASE X 


T HE 

PKO&L& 17\ 

/s ? 




/ A*} 

. da _ 



, 

Pe o» z x 


da 


91 *, o') = 

° ; 

* >c> 

Csic} 


<jCo,t>) = 

; 

O ^ o 

{ A*is 5 

f!Xm, 


finite. 


Cmf') 


r/f£//v(, the cAPtAcE. nt/tmeoe. m o (A,) 4 a/J> 

e #0 / f T/oM CU-) GWjtS ; 

eg? - P C ±01 _ » 4 ^ . a 

E (pUAr/o/J (is) tfjs A j>o*ur-tnM Of FoPM: 

gt%s) = c, erpCP.z) + Cz.&'fptf’**) 

W *£■££■ r—i 

/> 6%/a . ) - 1 ■ - yp&/4 ■+ /SfJ 

P A = (Pk/Z.) - ifi£f4 i- &s 

ApPOCA-TiOA/ OF- eqUATiGhJ (AO) t 7'Z/}A/Z,EogA\&b t & (!/£.£ '. 

C , ~c> 

APPC.{C.4'rto/j 6(C yMC Tj£/;/us/=oerh£A Sou/j£)Aet eo/j£>/r/o A/ CJU^) 

C /</££: 

C4 = i) 

~7tf£Z£FoZ£ } t#E 7-/C4A/sF£:£, E’o/JC T~(CAJ nP- me LAPLAC-E £>Ofl\AiN 

twpULse £&=,j>o*js& /s; 

%£ *>*) “ ex P [[£*£) ~i(i£/4) t &7J zj 
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rp& OOTL£r } TH-tS £-&I>Oc££ ~TC sr 

[(fe/z)- HZ p o 40 + ?cs J (va) 

7HE ££ LA-7/oA)*&ff> Po/L 7 fi£ P/iesT ftom&PT ABcuy- 1»£ 

omt, /// /s: 



PfEf oLsn/Ai 6 7~Af£. P/^fZ&e^/V T/Art^A} 4-/V£> Tj+£/t/(> T(+£ C/rtl/T 

VIZIRS’- 

Af. - X (VM) 

r 

T/te. Z-£ L A T/oAjSPtP Poe- TtHL ‘SBCvVt) /VO/rfE/Vr PSoUT TH£ 
flt£AV fs ; 



P£e^o^(7)/4/ c rPE X>/AP£L££?/JT'A-TiOfiJ } irAPIH6> TH£. V-tm/r, A Ad P^PXoOn /A}G- 
t&£ oofrrZAcTicP yiexps: 



C^sUo) 


s* zzr 


TAP P£oBc£rn vs' 

/ . <fy 4M 

Pe ■ ~ d6 

ij &&,€>)= c> 4 Z*o 

d ^ (z 0) &)/$ g = o ^ & pc. 


Cax) 


OO 

{ Ka) 


y to,*} = As - a *0 


CV-aI) 


TAktfijC* -TAZ £APv-Ac£ TPAMzPoZM oP Cyix) AaO APP^Vt*6- 
eqo* Pro// ( X) C/V£S; 


d l __ ctyc^s) 

ctz 


sPg^ 5 ; 


^? 2 - 


= o 


( Xs**^) 
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77 5OLVT10 A JO Cx^ti J fS OF JAP FoEpl: 

< 3 (*,s) - etpCP,*) + c z z) 

ujpzee 

P, = * Y(Pc/+) + P*s ' 

P# ~ (fete) - rfcPe/ 4 ) + fe-T 
fif > f > LHltl& TtfE. T£AtV ( , FO0/V} oF £ QOArfOfi/ fak) £>/U£S : 

o - c,{~ e fP [( p -s 4 Jofai-fe' 

* ^ ( # --]pAT^T) etffffy* - (x// ) 

A PP at fit T#Z TZfivsFotm op eQJAnoto t)UO)j (A*!), C /</£$: 

± ~ £,+£?_ (xv) 

£q>v 4 r/o/Vs Or/v) 4 m (k/) A££- two pquAt/o/us /a/ two uajmow/js) 

SoLVli 06 FO<L C, AV6 G z /UVJ> S<*>&%lTt7VT/AlQ /A/JCz T* E- eUPXJSLZ^nft 

FoC tyC^s) y/F*J>s : 
ft%si = gyp ( p ji) 

{(§ A ®yp [ ? “V| i+ ^ 5 j ~ ( 3 ? ~W + fc ) | 

?*x, T#f£. Becomes: 




eyp 


)exp[-6-Z 

»if+ **j 

] 4-i 

~+Pes jeyp[o 

■4M{ 

i 

i 

i j 

^p{ 



hW^\ 

!{ 




AS: 


KSLATfoMSAtP Fofi- rm F/ZS.T Aloft EA)T ASotjT 



— j(Uw<- 


5 -* C. 


7 

<*r J 


Pstopm/Au Tfts- j> fFp^pp/j npncjj #a/j> 


TtfE OP/<UAJ 


FFtZAUg TH£ Ctft/r y /££.&$• 
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//, r=T X + [cypt-Pe?*)- «Yfi{Pe-Pe*-y) Jp<2 WM*) 


T#e k&cat/omzh/p tot- t#£ seeow #iom£A/r A&oor 

T#E ShEA/J /s: 



P&efop yn/#c THE PtEFEPe/VTlAT/W) TAJzwl TA£ Ltryj/T ) A*Jt> 
?S£FO£ /V/A/C THE So&7¥AtTiC>H ! t £U>S '• 


= *1?* + e Y P t Pe - Pe [(4/re ) - 6 &>/Pe) - &PeT ] 

+ [«¥p C- J Pe ^/CPt'f -\etP (1 Pe - 2 Pe a 1 j/(Pe?~ 

. . . Cxyjii) 
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Appendix F Upper and. Lower Bounds for the Condition of a Matrix* 
t/PP£& AN<£> cotfeZ £01 wtS om TNG- doNbirio*/ oG A /n/*Fr£iX_ 


77/£ CO/UA / TioAJ of A A>4T£tt Mrrtf Z&9t>&cr 7© /NVe Zs/oA/ 

ANA 4 PAWMtAe- A/oErt] h Pef/A/SA AS' 

tcCA) = II All -Ilf II 

A t$ TN£ rtjATfn 8 £/NC CeJJ $./£>£ . 

II II is ritz jntre/y Noam se/M# oseJ>. 

AfOOJ c£T" 

Y == A • u>tt&££ *Y is> the. <lo v aPut‘ei> f^ve^s £ A . 


&££:/#£ TN£ /W/tTjeiX e.£SfPOAL~j £ ’ 

£_ = AT - r ^ UIH£E£ X VS THE xPEWTVM 

Tt/B&BfoeiB. ; 

i = £’ + ='s 


AW 3 >: 


ALSO: 


IIY 1 I & liril-ljr-f |\l 
a - 1 = Y-tT-^ir 1 - 


W) 


TVViS £r*PUE£, 

HA -11 ^ \ 1 Y 1 I • l\tc+|T 11 ^ 

mUL'TPOUKJA EQuAtvoV /X} A^p W 4 \l CotABlMU & 

ft Mill £ UMHUT & 

ftxt£H 

° c: 

Hft-ivp ^ vctt) - MHllHKriLT'U ^ 



